Quantum phase transitions in a pseudogap Anderson-Holstein model 
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We study a pseudogap Anderson-Holstein model of a magnetic impurity level that hybridizes 
with a conduction band whose density of states vanishes in power-law fashion at the Fermi energy, 
and couples, via its charge, to a nondispersive bosonic mode (e.g., an optical phonon). The model, 
which we treat using poor-man's scaling and the numerical renormalization group, exhibits quantum 
phase transitions of different types depending on the strength of the impurity-boson coupling. For 
weak impurity-boson coupling, the suppression of the density of states near the Fermi energy leads 
to quantum phase transitions between strong-coupling (Kondo) and local-moment phases. For 
sufficiently strong impurity-boson coupling, however, the bare repulsion between a pair of electrons 
in the impurity level becomes an effective attraction, leading to quantum phase transitions between 
strong-coupling (charge-Kondo) and local-charge phases. Even though the Hamiltonian exhibits 
different symmetries in the spin and charge sectors, the thermodynamic properties near the two 
types of quantum phase transition are closely related under spin-charge interchange. Moreover, the 
critical responses to a local magnetic field (for small impurity-boson coupling) and to an electric 
potential (for large impurity-boson coupling) are characterized by the same exponents, whose values 
place these quantum critical points in the universality class of the pseudogap Anderson model. One 
specific case of the pseudogap Anderson-Ifolstein model may be realized in a double-quantum-dot 
device, where the quantum phase transitions manifest themselves in the finite-temperature linear 
electrical conductance. 

PACS numbers: 71.10.Hf, 75.30.Hx, 72.10.Di, 64.60.ae 



I. INTRODUCTION 

Quantum phase transitions (QPTs) take place between 
competing ground states at the absolute zero of temper- 
ature (T = 0) upon variation of a nonthermal control 
parameteriii^ QPTs are thought to play a role in many 
important open problems in condensed-matter physics, 
including high-temperature superconductivity;^^— the 
phase diagram for magnetic heavy-fermion metals^ii and 
various types of metal-insulator transitiouj^i^ 

An interesting class of zero-temperature transitions is 
impurity or boundary QPTs at which only a subset of 
system degrees of freedom becomes critical.^- A well- 
studied example arises in the pseudogap Anderson im- 
purity modeliir— of an interacting impurity level hy- 
bridizing with a host density of states that vanishes in 
power-law fashion precisely at the Fermi energy — a prop- 
erty that can be realized in a number of systems in- 
cluding unconventional d-wave superconductors,— cer- 
tain semiconductor heterostructures,^- and a particu- 
lar double-quantum-dot setup»^"— The reduction of the 
density of states near the Fermi energy leads to QPTs 
between strong-coupling (Kondo-screened) and local- 
moment phaseSfU. At the transitions, the system exhibits 
a critical response to a local magnetic field applied only 
to the impurity sitei^"— 

Impurity quantum phase transitions have been 
predicted— ^^^1^^" — and possibly observed^^^ to arise 
in various nanodevices. While strong electron-electron 
interactions are an integral element of such nanode- 
vices, experiments on single- molecule transistors'^^ and 
quantum-dot cavities^ have also highlighted the impor- 
tance of electron-phonon interactions. The main aspects 



of the last two experiments appear to be captured by 
variants of the Anderson-Holstein model, which supple- 
ments the Anderson model^l for a magnetic impurity in 
a metallic host with a Holstein coupling'^* of the impu- 
rity charge to a local bosonic mode, usually assumed 
to represent an optical phonon. The model has been 
studied since the 1970s in connection with the mixed- 
valence problem;^"— the role of negative-[/ centers in 
superconductors )"^ '''"^^ and most recently, single-molecule 
deviceSi^"^ Various analytical approximations as well as 
nonperturbative numerical renormalization-group calcu- 
lations have shown that the Holstein coupling reduces the 
Coulomb repulsion between two electrons in the impurity 
level, even yielding effective electron-electron attraction 
for sufficiently strong impurity-boson coupling. Further- 
more, for the full Anderson-Holstein model with nonzero 
hybridization, as the impurity-boson coupling increases 
from zero, there can be a crossover from a conventional 
Kondo effect, involving conduction-band screening of the 
impurity spin degree of freedom, to a "charge Kondo ef- 
fect" in which it is the impurity "isospin" or deviation 
from half-filling that is quenched by the conduction band. 
However, the evolution between these limits is entirely 
smooth, and the model exhibits no QPT. 

This paper reports the results of study of a pseudogap 
Anderson-Holstein model which incorporates the struc- 
tured conduction-band density of states from the pseudo- 
gap Anderson model into the Anderson-Holstein model. 
The essential physics of the problem, revealed using a 
combination of poor-man's scaling and the numerical 
renormalization group (NRG), is shown to depend on 
the sign of the effective Coulomb interaction between two 
electrons in the impurity level, on the presence or absence 
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of particle-hole (charge-conjugation) symmetry and time- 
reversal symmetry, and on the value of the exponent r 
characterizing the variation p{£) cx jej'" of the density of 
states near the Fermi energy £ = 0. Even though the 
pseudogap Anderson-Holstein Hamiltonian has a lower 
symmetry than the pseudogap Anderson Hamiltonian, 
the universal properties of the former model, including 
the structure of the renormalization-group fixed points 
and the values of critical exponents describing properties 
in the vicinity of those fixed points, are identical to those 
of the latter model as generalized to allow for negative 
(attractive) as well as positive values of the local inter- 
action U between two electrons in the impurity leveli^ 
The pseudogap Anderson-Holstein model can therefore 
be regarded in part as providing a physically plausible 
route to accessing the negative-C/ regime of the pseudo- 
gap Anderson model. Anderson impurities with [/ < 
have recently attracted attention as a possible route to 
achieving enhanced thermoelectric power.-- 

The remainder of this paper is organized as follows: 
Section introduces the pseudogap Anderson-Holstein 
model, analyzes special cases in which the model reduces 
to problems that have been studied previously, outlines a 
perturbative scaling analysis of the full model, and sum- 
marizes the numerical renormalization-group approach 
used to provide nonperturbative solutions of the model. 
Section IHII presents results under conditions of particle- 
hole and time-reversal symmetry while Sec. IIVI addresses 
the general model with band exponent r between and 
1. Section|V]focuses on the specific case r = 2 relevant to 
a boson-coupled double-quantum-dot device. Section IVll 
summarizes the main results of the paper. The Appendix 
El contains details of the perturbative scaling analysis. 



II. MODEL HAMILTONIAN, PRELIMINARY 
ANALYSIS, AND SOLUTION METHOD 

A. Pseudogap Anderson-Holstein model 

In this work, we study the pseudogap Anderson- 
Holstein model described by the Hamiltonian 

a — -^imp^"-^band 4" -^boson^" -^imp-band ^ -^imp-boson; (1) 

where 



-^imp 
^^band 



- 5d [fid - 1) + \U{hd - if 



kcri 



-^boson 
-^imp-band 



k,o- 



Hi 



imp-boson 



;7^E(^k4ck. + H.c.), 
X{hd-l){b + b^). 



(2a) 
(2b) 

(2c) 
(2d) 

(2e) 



Here, da annihilates an electron of spin z component 
a = ±i (or a = f , J,) and energy Sd = Sd — \U < 



in the impurity level, fid = X^cr '^da (with fida = f^J-f^cr) is 
the total impurity occupancy, and i7 > is the Coulomb 
repulsion between two electrons in the impurity level 
Vk is the hybridization matrix element between the im- 
purity and a conduction-band state of energy Sk anni- 
hilated by fermionic operator Cko-, and A characterizes 
the Holstein coupling of the impurity occupancy to the 
displacement of a local vibrational mode of frequency 
ujQ. Nk is the number of unit cells in the host metal 
and, hence, the number of inequivalent k values. With- 
out loss of generality, we take Vk and A to be real and 
non-negative. For compactness of notation, we drop all 
factors of the reduced Planck constant h, Boltzmann's 
constant ks, the impurity magnetic moment QHb, and 
the electronic charge e. 

The conduction-band dispersion ek and the hybridiza- 
tion Vk affect the impurity degrees of freedom only 
through the hybridization functioit^ 



(3) 



To focus on the most interesting physics of the model, 
we assume a simplified form 



T{e)^T\e/D\'-e{D-\e\), 



(4) 



where Q{x) is the Heaviside function and we refer to the 
prefactor F as the hybridization width. In this notation, 
the case r = represents a conventional metallic hy- 
bridization function. This paper focuses on cases r > 
in which the hybridization function exhibits a power-law 
pseudogap around the Fermi energy. One way that such 
a hybridization function can arise is from a purely local 
hybridization matrix element Vk = V combined with a 
density of states (per unit cell per spin orientation) vary- 
ing as 

Pie) ^ iV-i J2 - ^k) = Pole/Dl-eiD - \s\), (5) 



in which case F = irpoV"^. However, the results presented 
in this paper apply equally to situations in which the k 
dependence of the hybridization contributes to the energy 
dependence of T{e). 

The assumption that r(e) exhibits a pure power-law 
dependence over the entire width of the conduction band 
is a convenient idealization. More realistic hybridization 
functions in which the power-law variation is restricted to 
a region around the Fermi energy exhibit the same qual- 
itative physics, with modification only of nonuniversal 
properties such as critical couplings and Kondo temper- 
atures. 

The properties of the Hamiltonian specified by Eqs. 
turn out to depend crucially on whether or not 
the system is invariant under the particle-hole transfor- 
mation Cj.^ — > cj^^, dg. —d^., b — >■ —b, which maps 
ffk — £k and 6d — > — <5d- For the symmetric hybridiza- 
tion function given in Eq. Q, the condition for particle- 
hole symmetry is (5^ = corresponding to Sd 



-\U. 
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B. Review of related models 

Before addressing the full pseudogap Anderson- 
Holstein model, it is useful to review two limiting cases 
that have been studied previously. 

1. Pseudogap Anderson model 

For coupling A = 0, the pseudogap Anderson- 
Holstein model reduces to the pseudogap Anderson 
modeli2r2ii2ii2^ plus free local bosons. In the conven- 
tional (r = 0) Anderson impurity model, the generic low- 
temperature limit is a strong-coupling regime in which 
the impurity level is effectively absorbed into the con- 
duction band.-^ In the pseudogapped (r > 0) variant 
of the model, the depression of the hybridization func- 
tion around the Fermi energy gives rise to a competing 
local-moment phase in which the impurity retains an un- 
screened spin degree of freedom all the way to absolute 
zero. The T = phase diagram of this model depends 
on the presence or absence of particle-hole symmetry^* 
and of time-reversal symmetryi^^ 

Behavior at particle-hole symmetry (Sd — 0): For any 
band exponent < r < ^, in zero magnetic field there is 
a continuous QPT at a critical coupling F = Fc(r, U, Sd = 
0) between the local-moment phase and a symmet- 
ric strong-coupling phase. In the local-moment phase 
(reached for < F < Fc), the impurity contributions^ to 
the entropy and to the static spin susceptibility approach 
the low-temperature limits S'imp = ln2 and Txs,imp — 
1/4, respectively, while conduction electrons at the Fermi 
energy experience an s-wave phase shift (5o(e = 0) = 0. In 
the symmetric strong-coupling phase (F > F^), the cor- 
responding properties are S'imp = 2rln2, Txs.imp = r/S, 
and So{0) = —(1 — r)(7r/2) sgne, all indicative of partial 
quenching of the impurity degrees of freedom. 

A magnetic field that couples to the band electrons 
moves the zero in the density of states of each spin species 
away from the Fermi level and washes out all pseudogap 
physics at energies below the Zeeman scalci^ More in- 
teresting is the breaking of time-reversal symmetry by a 
local magnetic field that couples only to the impurity de- 
gree of freedom and enters the Anderson model through 
an additional Hamiltonian term 

h 

Hh^ ^(ndt~ndi)- (6) 

The critical response to an infinitesimal h reveals that 
the transition between the local-moment and strong- 
coupling phases takes place at an interacting quantum 
critical point.-^-"^^ However, a finite value of h destabi- 
lizes both the local-moment phase^^ and the symmetric 
strong-coupling phasejS and destroys the QPT between 
the two.^^ For any h ^ 0, the ground state of the particle- 
hole-symmetric model (with > 0) is a fully-polarized 
local moment that is asymptotically decoupled from the 
conduction bandit 




FIG. 1: Schematic 5d-T phase diagrams of the pseudogap An- 
derson model [Eqs. (O-Q with A = 0] for band exponents (a) 
< r < i, (b) r > ^. Generically, the system falls into either 
a local-moment phase (LM) or one of two asymmetric strong- 
coupling phases (ASC±). However, there is also a symmetric 
strong-coupling phase (the line labeled SSC) that is reached 
only for < r < i under conditions of strict particle-hole 
symmetry {Sd ~ 0) and for sufficiently large hybridization 
widths r. 

For r > ^, the symmetric strong-coupling fixed point is 
unstable even in zero magnetic field ) ^'^d^ and a particle- 
hole-symmetric system lies in the local-moment phase for 
all values of F. 

Behavior away from particle-hole symmetry {dd 7^ 0): 
In zero magnetic field, the model remains in the local- 
moment phase described above for all l^^l < (i.e., 
-U < Ed < 0) and F < T,{r,U,Sd) = r^{r,U,-Sd). 
As shown schematically in Fig. [U the critical hybridiza- 
tion width Fc increases monotonically from zero as 
\Sd\ drops below ^U. For < r < i [Fig. Ufa)], 
Fc(r, U, 6d) smoothly approaches the symmetric critical 
value Fc(r,C/,0) as Sd ^0. For r > i [Fig. [^b)], 
Fc(r, U, Sd) instead diverges as Sd 0, consistent with 
the Sd = behavior discussed above. 

For Sd and F > Fc, the model lies in one of two 
asymmetric strong-coupling phases that share the low- 
temperature properties S'imp = and Txs^imp — 0. For 
Sd > 0, the Fermi-energy phase shift is <5o(0) — — Trsgne, 
while the ground-state charge (total fcrmion number 
measured from half- filling) is Q = —1. For Sd < 0, hy 
contrast, So{0) — +7rsgne and Q = +1. We label these 
two phases ASC_ and ASC+ according to the sign of Q. 

For r < r* ~ 3/8, the low-temperature physics on the 
phase boundary F — Fc(r, U,5d ^ 0) is identical to that 
at F = Fc(r, t/, 0), whereas for r > r* the properties 
are distinct J^i^ For r* < r < 1, the response to an 
infinitesimal local magnetic field shows that asymmetric 
transitions take place at two interacting quantum critical 
points (one for Sd > 0, the other for Sd < 0). For r > 
1, the QPTs are first-order— and can be interpreted as 
renormalized level crossings between the local-moment 
doublet and the ASC± singlet ground states. 

The asymmetric strong-coupling phase is stable over 
a range of local magnetic fields. However, at a critical 
value oi \h\ the system undergoes a level-crossing QPT 
into the same fully polarized phase as is found at particle- 
hole symmetry^ 
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Relationship to the pseudogap Kondo model: In cases 
where F ^ ^C^ — li^dl, the pseudogap Anderson model can 
be mapped via a SchriefFer- Wolff transformation^ onto 
the pseudogap Kondo modeL— The latter model exhibits 
QPTs entirely equivalent to those described above~4 
This allows us to identify critical exponents obtained pre- 
viously for the pseudogap Kondo model^S as the values 
that apply to the special case A = of the pseudogap 
Anderson-Holstein model. 

2. Anderson-Holstein model 

For r = 0, the pseudogap Anderson-Holstein model 
reduces to the Anderson-Holstein modeli^"— Insight into 
the physics of both models can be gained by performing 
a canonical transformation of the Lang-Firsov type^ to 
eliminate the Holstein coupling between the bosons and 
the impurity occupancy [Eq. (I2el) ]. The transformation^'^ 

H^e'^He-'^ with = — (n^ - - 6) (7) 

(jjQ 

maps Eq. ^ to 

— -^imp ^ -^band ^~ -^boson ^ -^imp-band; (8) 

in which boson and iJband remain as given in Eqs. ([2c)) 
and (|2b)) . respectively. Himp is identical to i/imp [Eq. 
([2a|) ] apart from the replacement of U by 

U^U-2ep, (9) 

where the polaron energy 

Bp = xyujo (10) 

represents an important energy scale in the problem. The 
invariance of 6d under the mapping implies a renormal- 
ization of the level energy from ej, — 5^ — to 

£d^£d + £p- (11) 
Finally, the impurity-band coupling term becomes 

5i„.p-ba„d = 5] 14 (5^4ck. + H.c) , (12) 

with 

^ ^ g-(A/c.„)(6t_6)^ (13) 

The canonical transformation in Eq. ([7]) maps the local 
boson mode to 

6 = e^6e-^-6-(A/c^o)("d-l), (14) 

effectively defining a different displaced-oscillator basis 
for each value of the impurity occupancy Ud, namely, the 
basis that minimizes the ground-state energy of iJimp + 



J^boson + ^^imp-boson- The elimination of the Holstein cou- 
pling is accompanied by two compensating changes to 
the Hamiltonian: a reduction in the magnitude — or even 
a change in the sign — of the interaction within the im- 
purity level, reflecting the fact that Eq. pc| lowers the 
energy of the empty and doubly occupied impurity con- 
figurations relative to single occupation; and incorpora- 
tion into the impurity-band term [Eq. (I12p ] of operators 
B and that cause each hybridization event to be ac- 
companied by the creation and absorption of a packet of 
bosons as the local mode adjusts to the change in the 
impurity occupancy n^. 

The analysis of Eq. ([S]) is trivial in the case F = 
of zero hybridization where the Fock space can be parti- 
tioned into subspaces of fixed impurity occupancy rid = 
0, 1, and 2, and the ground state within each sector cor- 
responds to the vacuum of the transformed boson mode. 
It can be seen from Eq. © that the effective on-site 
Coulomb interaction changes sign at A = Aq, where 

Ao = v/^oC//2 . (15) 

For weak bosonic couplings A < Aq, the effective interac- 
tion is repulsive, and for \Sd\ < \U the impurity ground 
state is a spin doublet with Ud = 1 and cr = For 
A > Ao, by contrast, the strong coupling to the bosonic 
mode yields an attractive effective on-site interaction and 
for \5d\ < —\U the two lowest-energy impurity states 
are spinless but have a charge (relative to half filling) 
(3 = nrf — 1 = ±1; these states are degenerate only under 
conditions of strict particle- hole symmetry {5d — 0). 

Various limiting behaviors of the full Anderson- 
Holstein model with F 7^ are understood4^i^ 

(i) If loq and A are both taken to infinity in such a way 
that Sp defined in Eq. (flUl) approaches a finite value, 
the model behaves just like the pure-fermionic An- 
derson model with U replaced by U while F and 
5d are unaffected by the bosonic coupling (implying 
that Ed — 5d — \U is replaced by Ed)- 

(ii) In the instantaneous or anti-adiabatic limit F ^ 
Wo < 00, the bosons adjust rapidly to any change in 
the impurity occupancy; for uJo,U ^ £j,, the physics 
essentially remains that of the Anderson model with 
U U, while for u!o,U <C £p, there is also a re- 
duction from F to Fog = F exp(— Cp/wo) in the hy- 
bridization width describing scattering between the 
Ud = and Ud — 2 sectors, refiecting the reduced 
overlap between the ground states in these two sec- 
tors. 

(iii) In the adiabatic limit F ^ ljq, by contrast, the 
bosons are unable to adjust on the typical time scale 
of hybridization events, and neither U nor F under- 
goes significant renormalization. 

(iv) In the physically most relevant regime F < wo < 
U,D, NRG calculationsi^i^ show that for F < 
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\5d\ ^ U, there is a smooth crossover from a conven- 
tional charge-sector Kondo effect for A ^ Aq (and 
thus r <C [/) to a charge-sector analog of the Kondo 
effect for A Ao (and T —U). The primary goal 
of the present work is to explain how this physics 
is modified by the presence of a pseudogap in the 
impurity hybridization function. 



such intermediate states leads to the scaling equations 



C. Poor-Man's Scaling 



As a preliminary step in the analysis of the pseudogap 
Anderson-Holstein Hamiltonian, we develop poor-man's 
scaling equations describing the evolution of model pa- 
rameters under progressive reduction of the conduction 
bandwidth. Haldane's scaling analysis^^ of the metallic 
(r = 0) Anderson model in the limit U ^ D has previ- 
ously been extended to the pseudogap case r > 0, both 
for infinite^'^ and finite^^ on-site interactions U. Here, the 
analysis is further generalized to treat the anti-adiabatic 
regime of the Anderson-Holstein model with both metal- 
lic and pseudogapped densities of states. 

Our analysis begins with the Lang-Firsov canonical 
transformation ([T]). In the anti-adiabatic regime, it is 
a good approximation to calculate all physical proper- 
ties in the vacuum state of the transformed boson mode 
defined in Eq. ([HI) . We therefore focus on many-body 
states representing the direct product of the bosonic vac- 
uum with the half-filled Fermi sea and with one of the 
four possible configurations of the impurity level. In the 
atomic limit F = 0, the energies of the states having im- 
purity occupancy = 0, 1, and 2 can be denoted Eq, 
Ei=Eo+ Sd, and E2 = Ei + ea + U = 2Ei - Eq + U, 
respectively. 

The poor-man's scaling procedure involves progressive 
reduction of the conduction-band halfwidth from D to 
D. At each infinitesimal step D —i' D + dD < D, the 
energies Eq, Ei, and E2, as well as the hybridization 
function F are adjusted to compensate for the elimination 
of virtual hybridization processes involving band states 
in the energy windows D + dD < e < D and —D<e< 
— {D + dD). An added complication in the Anderson- 
Holstein model is the presence of the operators B and 
in Eq. (IT^ . which allow virtual excitation of states 
having arbitrarily high boson occupation numbers fib — 
b^b . As detailed in the Appendix [K[ summation over all 



dU _ 2F 
'dD~~ 



1 



1 



E(D+ed) £{D-ed) 
1 1 



d£d ^ F 

db TT 



8{D-U~ed) £{D + U+ed) 
1 2 



1 



SiD-id) SiD+id) SiD+U+id) 



dr 

dD 



T 
r — 
D 



(16) 

(17) 
(18) 



for renormalized model parameters that take bare values 
tj ^ U [Eq. ©], Ed = Ed [Eq. (HB)], and f = F for 
D = D. In these equations, the energy scale 



£iE)=E/S[-,^ 



is defined in terms of a dimensionless function 



S'(a, a;) = a e-l^l^ 



1 x^' 



n\ a + n 



aF(a) 7* (a, —x) e' 



(19) 



-\x\ 



(20) 

where F(a) is the gamma function and 7* (a, x) is related 
to the lower incomplete gamma function)^ 

In the case A = where £{E) = E, Eqs. (II11)-([IH1) re- 
duce to the scaling equations for the pseudogap Anderson 
modeli^ The pseudogap in the hybridization function 
produces a strong downward rescaling of F [see Eq. (|18p ] 
that leads, via Eqs. and ([TT]) to weaker renormaliza- 
tion of U and id than would occur in a metallic (r = 0) 
host. For A > 0, one finds that \£{E)\ > \E\, so the 
bosonic coupling acts to further reduce (in magnitude) 
the right-hand sides of Eqs. ([TE)) and ([T7| . and produces 
still slower renormalization of U and id with decreasing 
D. It should be noted that neither the bosonic energy ujq 
nor the electron-boson coupling A is renormalized under 
the scaling procedure, and that the scaling equations re- 
spect particle-hole symmetry in that bare couplings sat- 
isfying Ed — —\U inevitably lead to rescaled couplings 
that satisfy id = —\U. 

Equation (|18p can readily be solved to give 



t(D) = [b/DY F. 



(21) 



For A > 0, it is not possible to integrate Eqs. ([T6| and 
(IT71) in closed form due to the presence of the nontrivial 
function £{E) on their right-hand sides. The equations 
have been derived only to lowest order in nondegenerate 
perturbation theory, and are therefore limited in validity 
to the range I [/ 1, \id\, F ^ D. Nonetheless, one may be 
able to obtain useful insight into the qualitative physics of 
the model through numerical integration of Eqs. (|16p and 
(|17p until one of the following conditions is met, implying 
entry into a low-energy regime governed by a simpler ef- 
fective model than the full pseudogap Anderson-Holstein 
model: 
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(i) If id, U + 2ed > D > r, the system should enter the 
empty-impurity region of the strong-couphng phase, 
where the impurity degree of freedom is frozen with 
an occupancy close to zero. 

(ii) If -(C/ + Ed), -{U + 2ed) > D > the system 
should enter the full-impurity region of the strong- 
coupling phase, where the impurity degree of free- 
dom is frozen with an occupancy close to two. 

(iii) If —id, U + id> D > T, the system is expected to 
enter an intermediate-energy local-moment regime 
in which the impurity states with rt^ ^ 1 are 
frozen out. As discussed further in Sec. IIII Bl one 
can perform a generalization of the SchriefFer- Wolff 
transformation^ to map the pseudogap Anderson- 
Holstein model to a pseudogap Kondo model with 
the density of states in Eq. ([S]). Depending on the 
value of the Kondo exchange coupling generated 
by the SchriefFer- Wolff transformation, the system 
may lie either in the strong-coupling phase of the 
pseudogap Kondo model (which should correspond 
to another region of the strong-coupling phase of 
the Anderson-Holstein model) or in a local-moment 
phase where the impurity retains a free two-fold spin 
degree of freedom down to absolute zero. 

(iv) l{id,-{U+id) > D >f, U+2id, the system should 
enter an intermediate-energy local-charge regime in 
which the = 1 impurity states become frozen 
out. A generalized Schrieffer- Wolff transformation 
can map the pseudogap Anderson-Holstein model to 
a pseudogap charge-Kondo model (see Sec. IIIIBI) . 
The system may lie either in the strong-coupling 
phase of the charge-Kondo model (yet another re- 
gion of the strong-coupling phase of the Anderson- 
Holstein model) or in a local-charge phase of both 
models (where the impurity retains a free two-fold 
charge degree of freedom down to absolute zero). 

(v) If f > i) > \id\ and/or f > D > \ij + id\, then the 
system should enter the mixed- valence region of the 
strong-coupling phase. 

Since each of the crossovers described above lies be- 
yond the range of validity of the scaling equations, the 
preceding analysis is only suggestive. In order to provide 
a definitive account of the pseudogap Anderson-Holstein 
model, it is necessary to obtain full, nonperturbative so- 
lutions, such as those provided by the numerical renor- 
malization group. However, we shall return to the scaling 
equations in Sec. lVBl to assist in the quantitative analysis 
of numerical results. 



D. Numerical solution method 

We have solved the model Eq. ([T]) using the numeri- 
cal renormalization-group (NRG) method^'^'^ as ex- 
tended to treat problems with an energy-dependent hy- 



bridization functioufi^iii and ones that involve local 
bosons.'*'' Briefly, the procedure involves three key steps: 
(i) Division of the full range of conduction-band energies 
—D < < D into a set of logarithmic intervals bounded 
by Em = iDA^™ for TO = 0,1,2,..., where A > 1 is 
the Wilson discretization parameter. The continuum of 
states within each interval is replaced by a single state of 
each spin a, namely, the linear combination of states ly- 
ing within the interval that couples to the impurity, (ii) 
Application of the Lanczos procedure to map the dis- 
cretized version of ffband onto a tight-binding form^ 

oo 

^band = E E(/-/n-l,. + H-C), (22) 

n— <y 

where /o^ oc X^k^kCk^, and {/„,^, = Sn.n' 5a,a' ■ 
The hopping parameters t„ (with = 0) contain all in- 
formation about the energy dependence of the hybridiza- 
tion function T{e). (iii) Iterative solution of the problem 
via diagonalization of a sequence of rescaled Hamiltoni- 
ans 

OS on -band oson Eg,o (23) 

and 

Hn = VAHn-1 + % E(/L/jv-i,a + H.c.) - Eg,n, 

(24) 

for iV = 1, 2, 3, . . ., where Eq.n is chosen so that the 
ground-state energy of Hn is zero. Hn can be inter- 
preted as describing a fermionic chain of length -I- 1 
sites with hopping coefficients that decay exponentially 
along the chain away from the end (site n = 0) to which 
the impurity and bosonic degrees of freedom couple. The 
solution of Hn captures the dominant physics at energies 
and temperatures of order DK~^/'^. 

The NRG procedure is iterated until the problem 
reaches a fixed point at which the spectrum of Hn and 
the matrix elements of all physical operators between the 
eigenstates are identical to those of Hn~2- (The eigen- 
solution of Hn differs from that of Hn-i even at a fixed 
point due to odd-even alternation effects.—) In addition 
to the conduction-band discretization, two further ap- 
proximations must be imposed. First, the number of 
states on the fermionic chain grows by a factor of 4 at 
each iteration, making it impractical to keep track of all 
the many-body states beyond the first few iterations. In- 
stead, one retains just the many-particle states of 
lowest energy after iteration N, creating a basis of di- 
mension 4iVs for iteration N + 1. Second, the presence 
of local bosons adds the further complication that the 
full Fock space is infinite-dimensional even at iteration 
N = 0, making it necessary to restrict the maximum 
number of bosons to some finite number A";,- 

The NRG calculations reported in the sections that 
follow took advantage of the conserved eigenvalues of the 
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total spin-z operator 

1 1 ^ 

s. = - {d\d^ - 4d,) + [jy^^ - fifj , (25) 

ri=0 

and the total "charge" operator 

N 

Q = n,~l + Y, (/,!t/nt + fljni - 1) (26) 

ji=0 

to reduce the Hamiltonian matrix to block-diagonal form, 
thereby reducing the labor of matrix diagonalization. In 
the absence of a magnetic field, the Hamiltonian 
commutes not only with Sz, but also with the total spin 
raising and lowering operators 

S+ = d\d, + J2fltfni and S^ = {S+y, (27) 

n 

the other two generators of SU(2) spin symmetry. By 
analogy, one can interpret 

71 

(28) 

as the generators of an SU(2) isospin symmetry. Since 
[-ffimp-boson, 7^ 0, Hn docs not exhibit full isospin ro- 
tation invariance, even though this symmetry turns out 
to be recovered in the asymptotic low-energy behavior at 
each of the important renormalization-group fixed points. 
In order to treat spin and charge degrees of freedom on 
equal footing, we elected not to exploit total spin con- 
servation in our NRG calculations. However, in the sec- 
tions that follow we identify NRG states by their quan- 
tum number S wherever appropriate. 

Throughout the remainder of this paper, all energies 
are expressed as multiples of the half-bandwidth D = 1. 
Results are reported for the representative case of a 
strongly correlated impurity level having U — 0.5 cou- 
pled to a local bosonic mode of frequency wq — O-l- 
The NRG calculations were performed using discretiza- 
tion parameter A = 2.5 or 3, allowing up to iVf, = 40 
bosons, values found to yield well-converged results for 
the model parameters considered. The number of re- 
tained many-body states was chosen sufficiently large to 
eliminate discernible truncation errors in each computed 
quantity; unless otherwise noted, this goal was attained 
using Ns — 500. 

III. RESULTS: PARTICLE-HOLE-SYMMETRIC 
MODEL WITH BAND EXPONENT < r < § 

As reviewed in Sec. IIIB| the particle-hole-symmetric 
pseudogap Anderson model with a band exponent < 
r < i has a QPT at F = T^r, U) between local- moment 
and symmetric strong-coupling phases. In this section we 
investigate the changes that arise from the Holstein cou- 
pling of the impurity charge to a local boson mode. For 



bosonic couplings A < Aq [see Eq. (ITSt ]. we find that the 
low-energy physics of the pseudogap Anderson-Holstein 
model is largely the same as for the pseudogap Anderson 
model with U replaced by an effective value Ucs [defined 
in Eq. (|33l) below] that differs from C7 introduced in Sec. 
111621 A QPT at F = F,i(r, (7, A < Ao) ~ F,(r, U,s) ex- 
hibits universal properties indistinguishable from those 
at the critical point of the pseudogap Anderson model. 
For stronger bosonic couplings A > Aq, there is instead 
a QPT at F = rc2{r,U,X > Aq) between the symmetric 
strong-coupling phase and a local-charge phase in which 
the impurity has a residual two- fold charge degree of free- 
dom. The critical exponents describing the local charge 
response at F = Fc2 are identical to those characterizing 
the local spin response at F = Fd. 

All numerical results presented in this section were ob- 
tained in a zero or infinitesimal magnetic field for a sym- 
metric impurity with Ed = ~\U = —0.25, for a bosonic 
frequency wq = 0.1, and for NRG discretization parame- 
ter A = 3. 



A. NRG spectrum and fixed points 

The first evidence for the existence of multiple phases 
of the symmetric pseudogap Anderson-Holstein model 
comes from the eigenspectrum of Hjsi . This spectrum can 
be used to identify stable and unstable renormalization- 
group fixed points of the model. 



1. Weak bosonic coupling 

FigurclDJa) shows— for r = 0.4, A = 0.05 < Aq ~ 0.158, 
and seven different values of F — the variation with even 
iteration number TV of the energy of the first excited mul- 
tiplet having quantum numbers S = 1, Q — 0. For small 
values of F, this energy En at first rises with increasing 
N, but eventually falls towards the value i?LM = ex- 
pected at the local-moment fixed point corresponding to 
effective model couplings F = A = and U = oo. At this 
fixed point, the impurity nd ~ 1 doublet asymptotically 
decouples from the tight-binding chain of length -I- 1, 
leaving a localized spin-i degree of freedom and low- lying 
many-body excitations characterized by a Fermi-energy 
s-wave phase shift (5o(e = 0) = 0, identical to that at 
the local-moment fixed point of the pseudogap Anderson 
model (see Sec.|IIB2j). 

For large F, instead rises monotonically to reach a 
limiting value i^ssc — 1.119 characteristic of the symmet- 
ric strong- coupling fixed point, corresponding to effective 
couplings F = c» and U — X = 0. Here, the impurity 
level forms a spin singlet with an electron on the end 
(n = 0) site of the tight-binding chain. The singlet for- 
mation "freezes out" the end site, leaving free-fermionic 
excitations on a chain of reduced length N, leading to 
a Fermi-energy phase shift Sq{Q) — —(1 — r)(7r/2) sgne. 
This is the same phase shift as is found at the symmetric 
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FIG. 2: (Color online) Particle-hole-symmetric pseudogap 
Anderson-Holstein model near its spin-sector critical point 
Cs: (a) NRG energy En vs even iteration number N of 
the first excited multiplet having quantum numbers 5* = 1, 
Q = 0, calculated for r = 0.4, U = -2£d = 0.5, cjq = 0.1, 
A — 0.05 < Ao — 0.158, and seven values of F — Fci (with 
Fci — 0.3166805) labeled on the plot, (b) Schematic phase 
diagram on the F-T plane for A < Ao, showing the T — 
transition between the local-moment (F < Fd) and symmet- 
ric strong-coupling (F > Fci) phases. Dashed lines mark the 
scale Ti* of the crossover from the intermediate-temperature 
quantum-critical regime to one or other of the stable phases, 
(c) Crossover scale Ti = DA'^^^^ vs |F - Fci| in the local- 
moment and symmetric strong-coupling phases, showing the 
power-law behavior described in Eq. (|38p. Here, is the in- 
terpolated value of at which En in (a) leaves its quantum- 
critical range by crossing one or other of the horizontal dashed 
lines. 

strong-coupling fixed point of the pseudogap Anderson 
modeliii^ 

The local-moment and symmetric strong-coupling 
fixed points describe the large- A^ (low-energy 
physics for all initial choices of the hybridization width 
except r = Fci ~ 0.3166805, in which special case En 
rapidly approaches Ec — 0.6258 and remains at that 
energy up to arbitrarily large N. This behavior can 
be associated with an unstable critical point Cs sepa- 
rating the local-moment and symmetric strong-coupling 
phases. (The subscript "S" indicates that Cs separates 
phases having different ground-state spin quantum num- 
bers.) The critical point corresponds to the pseudogap 
Anderson-Holstein model with A = and T/U equaling 
some r-dependcnt critical value. 

Whereas the critical coupling Fd is a nonuniversal 
function of all the other model parameters (r, U , wq, and 
A), the low-energy NRG spectra at the local-moment, 
symmetric strong-coupling, and Cs fixed points depend 
only on the band exponent r and the NRG discretiza- 
tion parameter A. For given r and A, each spectrum is 
found to be identical to that at the corresponding fixed 
point of the particle-hole-symmetric pseudogap Ander- 
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FIG. 3: (Color online) Particle-hole-symmetric pseudogap 
Anderson-Holstein model near its charge-sector critical point 
Cc: (a) NRG energy En vs even iteration number N of the 
first excited state having quantum numbers S = Q = 0, calcu- 
lated for r = 0.4, U = -2ed = 0.5, ljo = 0.1, A = 0.2 > Ao ~ 
0.158, and seven values of F - Fc2 (with Fc2 ^ 0.6878956) 
labeled on the plot, (b) Schematic F-T phase diagram for 
A > Ao, showing the T — transition between the local- 
charge and symmetric strong-coupling phases and the scale 
T2 of the crossover from the quantum-critical regime to a sta- 
ble phase, (c) Crossover scale Tj* = _DA"^2 /2 |p _ p^^l in 
the local-charge and symmetric strong-coupling phases, show- 
ing the power-law behavior described in Eq. (I39|l . Here, A| 
is the interpolated value of N at which En in (a) leaves its 
quantum-critical range by crossing one or other of the hori- 
zontal dashed lines. 



son model. Not only can the spectrum be interpreted as 
arising from an effective boson coupling A = 0, but it ex- 
hibits the SU(2) isospin symmetry that is broken in the 
full pseudogap Anderson-Holstein model. 



2. Strong bosonic cow^ 



Figure |3l[a) plots the energy at even iterations of 
the first NRG excited state having quantum numbers 
5 = Q = 0, for r = 0.4, A = 0.2 > An ~ 0.158, and 
seven different F values. For T > Vc2 ^ 0.6878956, 
En eventually flows to the value -Esse — 1-119 iden- 
tified in the weak-bosonic-coupling regime, and exami- 
nation of the full NRG spectrum conflrms that the low- 
temperature behavior is governed by the same symmetric 
strong-coupling fixed point. 

For r < Fc2, En flows to zero, the value found at 
the local-moment fixed point. In fact, all the fixed- 
point many-body states obtained for A > Aq turn out 
to have the same energies as states at the local-moment 
fixed point. However, the quantum numbers of states in 
the A > Ao spectrum and the local-moment spectra are 
not identical, but rather are related by the interchanges 
S ^ I and Sz Iz- We therefore associate the A > Ao 
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FIG. 4: (Color online) Phase boundaries of the particle- 
hole-symmetric pseudogap Anderson-Holstein model in zero 
magnetic field: Variation with bosonic coupling A of the 
critical hybridization widths Fci and Vc2- Data for U = 
— 2Ed = 0.5, LtjQ — 0.1, and band exponents r — 0.1, 0.2, 
0.3, and 0.4 are scaled and offset for clarity: the quanti- 
ties plotted are Fci (A) /Fci (0) -I- lOr — 1 (empty symbols) and 
Fc2(A)/20Fci(0) -I- lOr - 1 (filled symbols). Solid lines show 
the prediction of Eq. ((34)), while dashed lines plot the form 
^(r)A^'^"''' suggested by Eq. ((37]) with a prefactor A{r) de- 
termined by fitting values of Fc2 for 0.3 < A < 0.4. 



spectrum with a local- charge fixed point, corresponding 
to r = A = and U = — oo, at which the impurity has 
a residual isospin-i degree of freedom. Like its local- 
moment counterpart, this fixed point exhibits a phase 
shift 6o{0) = 0. 

For F — Fc2, rapidly approaches and remains at 
the same critical value Ec as found for A < Aq and 
r = Fci. Once again, however, the many-body spectrum 
is related to that at the corresponding weak-bosonic- 
coupling fixed point by interchange of spin and isospin 
quantum numbers, leading to the interpretation of this 
fixed point as a charge analog Cc of the critical point of 
the particle-hole-symmetric pseudogap Anderson model. 



bridization width F^ of the corresponding pseudogap An- 
derson model. The lines represent analytical expressions 
for the phase boundaries that will be explained in the 
remainder of this subsection. 

In the particle-hole-symmetric pseudogap Ander- 
son model, the critical hybridization width Tc{r,U) 
can be established by performing a SchriefFer- Wolff 
transformation^^ that maps the problem to a pseudogap 
Kondo model with a Kondo exchange coupling 



8T fU 



(29) 



Here, the Kondo coupling in a conventional metal (r = 0) 
is multiplied by a factor {U/2DY that accounts for the 
irrelevance of the hybridization width under poor-man's 
scaling [see Eq. (fT8|) ] while neglecting the much weaker 
renormalization of the on-site interaction [Eq. (ITB))]. The 
critical coupling Jc of the particle-hole-symmetric pseu- 
dogap Kondo model satisfies poJc — jc{r), where jdr) — 
r for r < ^ (Ref. ITT| ) and jc{r) ^ oo for r i (Refs. 
113, Hi, andtig). Combining this condition with Eq. 
yields 



' \2D 



(30) 



It is important to note that an equivalent expression 
has been derived within the local-moment approach to 
the pseudogap Anderson model without reference to a 
SchriefFer- Wolff transformation [see Eq. (6.10b) of Ref. 
Il5j and has been verified via NRG calculations As 
such, Eq. (|30l) with a suitably chosen value of jdr) is 
applicable even for r approaching i where charge fluctu- 
ations for F ~ Fc are too strong to allow mapping to a 
Kondo model. We now consider how Eq. (pO| should be 
modified to describe the phase boundaries of the pseudo- 
gap Anderson-Holstein model. 



B. Phase boundaries 



Weak bosonic coupling 



Figure S] shows phase boundaries for the symmetric 
pseudogap Anderson-Holstein model, as established for 
U = —2ed — 0.5, and ujq — 0.1 by examination of the 
NRG spectrum. In the atomic limit F = we find a 
level-crossing transition between the local-moment and 
local-charge phases at A = Aq ~ 0.15812(1), a value 
in excellent agreement with the prediction of Eq. dTSl) . 
(Throughout this paper, a digit in parentheses following 
a number indicates the estimated nonsystematic error in 
the last digit of the number.) For each of four values 
of the band exponent r, the figure plots the critical hy- 
bridization widths Fci (open symbols, for < A < Aq) 
and (filled symbols, for A > Aq) normalized by the 
A = value of Fd , which coincides with the critical hy- 



For A ^ Aq (and hence C7 > 0), it has been shownSS 
that a generalized Schrieffer- Wolff transformation maps 
the particle-hole-symmetric Anderson-Holstein model to 
a Kondo model with a dimensionless exchange coupling 



PoJ 



4F 



y. 1 (AM)^»^ 



(31) 



representing a sum over virtual transitions of the impu- 
rity from occupation = 1 to = or 2, accompa- 
nied by excitation of different numbers «{, = 0, 1, . . . of 
bosonic quanta. To facilitate comparison with the cor- 
responding expression p^J = ST/ttU for the symmetric 
Anderson model without bosons, one can use Eq. (I^IH) to 
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recast Eq. (1311) in the form 



8r 
sr 



\2ujQ loq J 



(32) 



where, both for C/ > (as is the case here) and for C/ < 0, 



C/cff = 2f(|[7|/2) 



(33) 



with £{E) as defined in Eq. Equation sug- 

gests that C/off plays the role of an effective Coufomb re- 
pulsion in the low-energy many-body physics of the full 
Anderson-Holstein model, distinct from the quantity U 
[Eq. 0] that emerges from considering just the atomic 
limit r = 0. Like f/, C/off passes through zero at A = Aq. 
For fixed U , the ratio Ues/U evolves smoothly from 1 for 
A = to e^/2^^° (~ 12 in the case U/2ujo ^ 2.5 used in 
our calculations) for A = Ao to 2 for A — > oo. 

Extension of the analysis of Ref. i50i to the case of 
a pseudogap density of states leads to the conclusion 
that the critical hybridization width separating the local- 
moment and symmetric strong-coupling phases should 
satisfy 



state of the displaced harmonic oscillator that minimizes 
the electron-boson interaction in the sector n^; = and 
the corresponding ground state for Ud = 2 (see Sec. 

EH. 

A poor-man's scaling analysis of the anisotropic pseu- 
dogap Kondo model^ indicates that for S> |J_l|, 
the phase boundary is defined by a condition poJz — 
r\n{2Jz/ J±). Applying this condition to the pseudogap 
Anderson-Holstein model, carrying over Eq. p6p from 
the case r = 0, and assuming [by analogy with Eq. p9p ] 
that poJz oc \Ues/D\^~^ , yields 



rc2(A) 



C/cff 



D 



(37) 



where we have used UcS — 2U ~ —iSp — — 4A^/wo for 
A » Aq. The validity of Eq. (|37)) is questionable because 
the critical hybridization widths it demands are too large 
to justify mapping to a Kondo model. Nonetheless, the 
NRG data for each value of r plotted in Fig. S] follow a 
A2(2-r) dependence (dashed lines) quite closely for 0.2 < 
A < 0.4. (This power law must break down closer to the 
level crossing between the local-moment and local-charge 
phases because necessarily vanishes at A = Aq.) 



rc=4j.(r)i? — 



(34) 



The solid lines plotted in Fig. |4] show the boundaries 
predicted by Eq. p4p with numerical evaluation of C/off- 
The agreement with the NRG data points is excellent for 
all four values of r, and for A extending from zero almost 
all the way to Aq. 



2. Strong bosonic coupling 

Cornaglia et al. have demonstrated^^^ that the 
Anderson-Holstein model with A 3> Aq maps to a 
charge analog of the Kondo model in which the impu- 
rity isospin degree of freedom [the impurity (d-electron) 
parts of the operators defined in Eqs. ([28ll ] is screened 
by its conduction-band counterpart. The impurity-band 
isospin exchange is anisotropic, with a longitudinal cou- 
pling pqJz = 8r/7r|C/off| and a transverse coupling 



PoJi 



sr 



7r|C/| V^i^o' Wo 



-'"-/"VoJ., (35) 



where U^a and S{a,x) are defined in Eqs. (l33t and (l20t . 
respectively. Closer investigation shows that an approx- 
imation that is equivalent for large A but also remains 
valid much closer to Ao is 



-|;7eft|/2Wo 7- 



(36) 



The strong suppression of "charge-flip" scattering arises 
from the exponentially small overlap between the ground 



C. Crossover scales 

Aside from allowing the identification of 
renormalization-group fixed points and phase bound- 
aries, the eigenspectrum of can also be used to define 
temperature scales characterizing crossovers between the 
domains of infiuence of different fixed points. We focus 
on the smallest such scale, which describes the approach 
to one of the stable fixed points of the problem. 



1. Weak bosonic coupling 

With decreasing [F — Fci|, Em in Fig. ^a) remains 
close to its critical value Ec over an increasing num- 
ber of iterations before heading either to i?LM or to 
^-ssc- To quantify this effect, it is useful to de- 
fine threshold energy values E± where i^LM < < 
Ec < E+ < Esse- The passage of En below E^ 
(above iJ+) at some Ni — determined by interpolation 
of the NRG data at even integer values of N — can be 
taken to mark the crossover around temperature T* ~ 
£)A~^i /2 from an intermediate-temperature quantum- 
critical regime dominated by the unstable critical point 
Cs to a low-temperature regime controlled by the stable 
local- moment (symmetric strong-coupling) fixed point. 
This crossover scale is expected to vanish for F — >■ Fd , as 
shown schematically in Fig.^b). Figure^c) plots values 
of T* determined by the criterion E^ = 0.3, E^ = 0.8. 
These data are consistent with the relation 



T* cx |F-F,i| 



r^Fei, 



(38) 
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TABLE I: Critical exponents at the critical point Cs of the 
particle-hole-symmetric pseudogap Anderson-Holstein model, 
evaluated for three band exponents r. The critical exponents 
are defined in Eqs. (|38|) and (I42p . A number in parentheses 
indicates the estimated random error in the last digit of each 
exponent. 



r 




/3i 


l/<5i 


Xl 
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0.2 


6.22(1) 


0.15(1) 


0.02630(2) 


0.9488(2) 


5.85(6) 


0.3 


5.14(1) 


0.34(1) 


0.07364(1) 


0.8629(3) 


4.41(3) 


0.4 


5.84(1) 


0.90(1) 


0.1845(1) 


0.6885(2) 


3.95(5) 



where vi is the correlation-length exponent at the quan- 
tum critical point. The numerical value of vi is indepen- 
dent of the precise choice of the thresholds E±. What 
is more, different combinations of the model parameters 
r, C/, clIq, and A result in different critical couplings Fd, 
but ui depends only on the band exponent r. Values for 
three representative cases are listed in Table HI 

2. Strong bosonic coupling 

The passage of E^^ outside a range < En < £'+ can 
also be used to define a crossover scale near its charge- 
sector critical point. This scale Tj* is expected to vanish 
at the critical point according to 

T*oc|r-r,2r^ as r^r,2, (39) 

a behavior that is sketched qualitatively in Fig. ^h) and 
is confirmed quantitatively in Fig.[3l[c). For all the values 
of r and A that we have studied, the numerical values of 
i>i and 1^2 coincide to within our estimated errors. 

D. Impurity thermodynamic properties 

This section addresses the variation with temperature 
of the impurity contributions^^ to the static spin and 
charge susceptibilities and to the entropy. With the con- 
ventional definitions Txs.imp — {Sf) and Txc,imp = {Q'^), 
a symmetric impurity level isolated from the conduc- 
tion band (F = 0) has TxsSmp = 1/4 and Txc, imp = 
for [/off' > T, but Txcimp 1 and Txs,imp = for 
UcS ^ ~T, with [/off' as defined in Eq. ([55)) . Due to the 
factor of 4 difference between the local-moment spin sus- 
ceptibility and the charge susceptibility of a local charge 
doublet, it is most appropriate to compare Txs;nxi-p with 
\TXc;im.p- During the NRG calculation of these thermo- 
dynamic properties, A^s = 3 000 states were retained after 
each iteration. 



1. Weak bosonic coupling 

Figure [5] plots the temperature dependence of Txs, imp, 
\Txc,im.-p: and S'imp for r = 0.4, A = 0.05, and seven 



r = 0.4, X = 0.05 




FIG. 5: (Color online) Thermodynamic properties of the 
particle-hole-symmetric pseudogap Anderson-Holstein model 
near its spin-sector critical point Cs: Temperature depen- 
dence of the impurity contribution to (a) the static spin sus- 
ceptibility Xs.imp multiplied by temperature, (b) the static 
charge susceptibility Xc.imp multiplied by temperature, and 
(c) the entropy Simp, for r = 0.4, U = —2ed = 0.5, cjo = 0.1, 
A = 0.05 < Ao — 0.158, and the seven values of F — Fd 
labeled in the legend. Filled (open) symbols connected by 
guiding lines represent data in the local-moment (symmet- 
ric strong-coupling) phase, while thick lines without symbols 
show the critical properties at Cs. Ns — 3 000 states were 
retained after each NRG iteration. 



values of F straddling Fd. At high temperatures T ^ 
max([/cff, F), the properties lie close to those of the free- 
orbital fixed point (Txs.imp = iTxc,imp = 1/8 and 
•S'imp = In 4), irrespective of the specific value of F. 
However, the T — >■ behaviors directly reflect the ex- 
istence of a QPT at F = Fd. In the local- moment 
phase (F < Fd), the residual impurity spin doublet 
is revealed in the limiting behaviors Txs.imp = 1/4, 
■j^^Xc.imp = 0, and S'imp = In 2. In the symmetric 
strong-coupling phase (F > Fd), the impurity degrees 
of freedom are quenched to the maximum extent possi- 
ble given the power-law hybridization function, yield- 
ing Txs.imp = jTxcAmp = r/8 and Si^p = 2rln2. Ex- 
actly at F = Fd, the low-temperature properties are 
distinct from those in either phase: Txs, imp — 0.1348, 
i^Xcimp ^ 0.0158, and ^i^p ^ 0.694 ~ In 2. These val- 
ues vary with the band exponent r, but are independent 
of other model parameters such as U, ujq, and A, so they 
can be regarded as characterizing the critical point Cs- 
For all the r values that we have examined, the critical 
properties coincide with those at the corresponding criti- 
cal point of the pseudogap Kondo or Anderson models.— 

When F deviates slightly from Fd, the thermody- 
namic properties follow their critical behaviors at high 
temperatures, but cross over for T < Tj* to approach 
the values characterizing the local-moment or symmet- 
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FIG. 6: (Color online) Thermodynamic properties of the 
particle-hole-symmetric pseudogap Anderson-Holstein model 
near its charge-sector critical point Cc: Temperature depen- 
dence of the impurity contribution to (a) the static spin sus- 
ceptibility Xs.imp multiplied by temperature, (b) the static 
charge susceptibility T^c, imp multiplied by temperature, and 
(c) the entropy Simp, for r — 0.4, U = —2ed = 0-5, cjq = 0.1, 
A = 0.2 > Ao — 0.158, and the seven values of F — Fc2 labeled 
in the legend. Filled (open) symbols connected by guiding 
lines represent data in the local-charge (symmetric strong- 
coupling) phase, while thick lines without symbols show the 
critical properties at Cc. Ns — 3 000 states were retained 
after each NRG iteration. 



ity at Cs mirrors that of the charge (spin) susceptibility 
at Cc, while the entropy behaves in the same manner 
at both critical points. These observations are consis- 
tent with the equivalence of the NRG spectra at the two 
fixed points under interchange of spin and charge quan- 
tum numbers (see Sec. IIII 



E. Local response and universality class 

In order to investigate in greater detail the proper- 
ties of the spin and charge critical points (Cs and Cc 
in Fig. nop , it is necessary to identify an appropriate 
order parameter for each QPT. The symmetric strong- 
coupling and local-moment phases can be distinguished 
by their values (0 and ^, respectively) of the magni- 
tude KS'z)! of the total spin in a vanishingly small mag- 
netic field applied along the z direction. Similarly, the 
magnitude \{Q)\ of the total charge in the presence of 
an infinitesimal electric potential takes the value in 
the symmetric strong-coupling phase and 1 in the local- 
charge phase. However, the fact that Sz and Q are 
conserved quantities — i.e., that the pseudogap Anderson- 
Holstein Hamiltonian commutes with Sz and Q defined 
in Eqs. (pS)) and I^E^, respectively — prevents these candi- 
date order parameters from exhibiting nontrivial critical 
exponents. Instead, we must look to the impurity 
response to local fields in order to probe the quantum 
critical behavior. 



ric strong-coupling phase. The crossover temperature T* 
coincides up to a constant multiplicative factor with that 
extracted from the NRG spectrum (as described in Sec. 
IIII B|l and its variation with |r — Fdl yields, via Eq. (|38l) . 
a correlation-length exponent i^i (r) in agreement with the 
values listed in Table HI 



2. Strong bosonic coupling 

Figure [HI plots the temperature dependence of Txs, imp, 
■jT^Xc.imp, and S'lmp for r = 0.4, A — 0.2, and various 
r straddling the critical value Tc2- Like in the case 
of weak bosonic coupling, the T — >■ behaviors distin- 
guish the two stable phases: the properties TxsAm-p = 0, 
jTxcSmp — 1/4, and S'lmp = In 2 in the local-charge 
phase contrast with Tx^^imp = \Txc,irap — f/?> and 
S'lmp = 2r In 2 in the symmetric strong-coupling phase. 
Exactly at F = Fc2, Txs.imp ^ 0.0158, |:Txc,imp ^ 
0.1348, and S'lmp — 0.694, values that can be taken to 
characterize the critical point Cc. From the thermody- 
namics near Cc, one can extract a crossover scale Tj* 
that gives [via Eq. (p9| ] a correlation-length exponent V2 
identical to that determined from the NRG spectrum. 

Figures [5] and [S] illustrate the general property that the 
temperature dependence of the spin (charge) susceptibil- 



1. Weak bosonic coupling 

In the pseudogap Kondo and Anderson models, the 
critical properties manifest themselves^^ through the re- 
sponse to a local magnetic field h that couples only to 
the impurity spin as specified in Eq. ([6]). The order pa- 
rameter for the pseudogap QPT is the limiting value as 
/i — >■ of the local moment 

Mioc = (i(ndt-"di)), (40) 

and the order-parameter susceptibility is the static local 
spin susceptibility 

Xs,\oc = - lim — — . (41) 

Based on the similarities noted above between the 
pseudogap- Anderson critical point and the Cs critical 
point of the pseudogap Anderson-Holstein model (i.e., 
the properties of the phases on either side of each transi- 
tion, the NRG spectrum at the transition, and the value 
of the order-parameter exponent), we expect that the two 
QPTs also to share the same order parameter. Accord- 
ingly, the behaviors of M\ac and Xs,ioc in the vicinity of 
the critical hybridization width F = Fd should be de- 
scribed by critical exponents 71, (5i, and xi defined 
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as follows: 

Mioc(r <Tci;h^O,T = 0)oc (Td - rf' , (42a) 

Xs,ioc(r > r,i;T = 0) cx (r - r^.y^ (42b) 
Mioc(ft;r-r,i,r = o)cx|/i|i/*S (42c) 

Xs,ioc(T;r = r,i) cxT-"\ (42d) 

The preceding expectations are proved correct by NRG 
calculations, as demonstrated in Fig. [7] for r = 0.4 and 
A = 0.05, the case treated in Fig. [5] The critical ex- 
ponents extracted as best-fit slopes of log-log plots are 
listed in Table U for three values of the band exponent 
r. The values of individual critical exponents vary with 
r, but are independent of other Hamiltonian parameters 
{U, Wo, and A) and are well converged with respect to 
the NRG parameters (A, Ng, and iVf,). To within their 
estimated accuracy, the critical exponents for a given r 
obey the hyperscaling relations 

1 -\- Xi 

Si = , 2/?i = z/(l - xi), 7i=z/ixi, (43) 

1 — a;i 

which are consistent with the scaling ansatz 

for the nonanalytic part of the free energy at an interact- 
ing critical pointi^ 

2. Strong bosonic coupling 



r = 0.4, ;l = 0.05 




FIG. 7: (Color online) Static local spin response of the 
particle-hole-symmetric pseudogap Anderson-Holstein model 
near its spin-sector critical point Cs. Circles are NRG data 
for r = 0.4, U = -2ed = 0.5, ujq = 0.1, and A = 0.05, 
at or near the critical hybridization width Fci — 0.3166805. 
Straight lines represent power-law fits, (a) Static local spin 
susceptibility Xs,ioc(r = 0) vs F— Fci in the symmetric strong- 
coupling phase, (b) Local magnetization Miodh ^ 0, T = 0) 
vs Pel — r in the local-moment phase. Inset: Continuous van- 
ishing of Mioc{h — >■ 0, T = 0) as F approaches Fci from below, 
(c) Local magnetization Afioc(r = Tc\,T — 0) vs local mag- 
netic field h. (d) Static local spin susceptibility Xs,ioc(r = Fc) 
vs temperature T. 



We have seen above that the NRG spectrum and low- 
temperature thermodynamics at the Cc fixed point are 
related to those at the Cs fixed point by interchange of 
spin and charge degrees of freedom. One therefore ex- 
pects to be able to probe the critical properties via the 
response to a local electric potential (p that enters the 
model through an additional Hamiltonian term 

H^ = c^{hd-l). (45) 

Comparison with Eq. pap shows that cj) is equivalent to 
a shift in Sd (or e^). The order parameter should be the 
(/> — limiting value of the local charge 

Qioc - {fid - 1), (46) 

and the order-parameter susceptibility should be the 
static local charge susceptibility 

Xc.ioc = - lim (47) 

In the vicinity of the critical point F = Fc2, one expects 
the following critical behaviors: 

Qioc(r < Fe2; ^ 0, T = 0) cx (Fe2 - F)''^ , (48a) 

Xc,ioc(r>F,2;T = 0)« (F-F,2)"^% (48b) 

Qioc(0;F F,2,T = 0) (X (48c) 

Xc,ioc(r;F = Fe2)cxr-=^^ (48d) 



These expectations are borne out by the NRG results, as 
iUustrated in Fig. |5] for the case r = 0.4, A = 0.2 treated 
in Fig. 13 



3. Comparison between weak and strong bosonic coupling 

Figure Hl^a) superimposes the variation with F of the 
order parameter in the vicinity of the Gs and Cc critical 
points for two representative band exponents, r = 0.2 
and 0.4. The equality of the slopes of the log-log plots 
at the spin- and charge-sector QPTs shows that /?i = /32- 
Similarly, Fig. IH^b) shows that the temperature varia- 
tion of the order-parameter susceptibilities is consistent 
with xi — X2- Indeed, for each value of r that we have 
examined, we find that all critical exponents at Gc are 
indistinguishable (within our estimated errors) from the 
corresponding exponents at Cs and at the critical point 
of t he p seudogap Kondo model (as given in Table I of 
Ref . 1231 ) . This leads us to conclude that all three critical 
points lie in the same universality class. 



F. Renormalization-group flows 

The essential physics of the particle-hole-symmetric 
pseudogap Anderson-Holstein model can be summarized 
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r = 0.4, X = 0.2 




FIG. 8: (Color online) Static local charge response of the 
particle-hole-symmetric pseudogap Anderson-Holstein model 
near its charge-sector critical point Cc- Circles are NRG data 
for r = 0.4, U = —2ed = 0.5, wo = 0.1, and A — 0.2, at or near 
the critical hybridization width Tc2 — 0.6878956. Straight 
lines represent power-law fits, (a) Static local charge suscep- 
tibility Xc,ioc(r = 0) vs r — Tc2 in the symmetric strong- 
coupling phase, (b) Local charge Qioc(0 — >■ 0,r = 0) vs 
Tc2 — r in the local-charge phase. Inset: Continuous vanish- 
ing of (5ioc(0 — >■ 0,T = 0) as r approaches Tc2 from below, 
(c) Local charge Qioc(r = rc2,T = 0) vs local electric poten- 
tial (j). (d) Static local charge susceptibility Xc,ioc(r = rc2) vs 
temperature T. 



oM^jh^o,T=o) ox,,,jr = rei) 

□ Q,j(^^o,r=o) <x,ioc(r = rc2) 




FIG. 9: (Color online) Comparison of the static local re- 
sponses of the particle-hole-symmetric pseudogap Anderson- 
Holstein model near its critical points Cs and Cc: (a) De- 
pendence of the local magnetization M\oc{h — > 0, T = 0) on 
Fci — F and of the local charge Qioc{4> — >■ 0, T = 0) on Fc2 — F, 
for U = —2ed = 0.5, — 0.1, A = 0.05 (magnetic response) 
or 0.2 (charge response), and band exponents r — 0.2 and 
r = 0.4. (b) Static local spin susceptibility Xs,ioc(r; F = Fci) 
and static local charge susceptibility Xc,ioc(r;F = Fc2) vs 
temperature T. All parameters other than F take the same 
values as in (a). Straight lines represent power-law fits. 



in the schematic renormalization-group flow diagram 
shown in Fig. [lOl which applies to all band exponents 
in the range < r < i . Arrows indicate the evolution of 
the effective Coulomb interaction UoS and the hybridiza- 
tion width F with increasing NRG iteration number 
N, i.e., under progressive reduction of the temperature 
T ~ DA~^^^. The high-temperature limit of the model 
is governed by the free-orbital (FO) fixed point, corre- 
sponding to bare model parameters J7 = F = A = and a 
Fermi- level phase shift i5o(0) = 0. Dashed lines mark the 
separatrices between the basins of attraction of the local- 
moment (LM), local-charge (LC), and symmetric strong- 
coupHng (SSC) fixed points described above. Flow along 
each separatrix is from the free-orbital fixed point to- 
wards one or other of two quantum critical points — either 
the conventional spin-sector critical point Cs reached for 
Ucff > 0, or its charge analog Cc reached for [/off < 0. 

A renormalization-group fixed-point structure equiv- 
alent to that described in the preceding paragraph has 
been presented previousljs^^ for the pseudogap Ander- 
son model under the assumption that the bare on-site 
Coulomb interaction U may be taken to be positive or 
negative. Indeed, many of the universal properties of 
the pseudogap Anderson-Holstein model presented in this 
section — particularly ones associated with the quantum 
critical points Cs and Cc — reproduce those of this ex- 
tended pseudogap Anderson model. 



However, we emphasize that the particle- hole- 
symmetric Anderson and Anderson-Holstein models have 
different symmetries and are therefore not trivially re- 
lated to one another. The pseudogap Anderson Hamilto- 
nian exhibits exact SU(2) spin and isospin (charge) sym- 
metries, and all physical properties at a point {U,T) in 
the diagram analogous to Fig. [TU] [see Fig. 1(b) of Ref. 
[25| map exactly to the properties at {—U, F) under the 
interchange of spin and charge degrees of freedom. No 
such mapping holds in the pseudogap Anderson-Holstein 
model, where the Hamiltonian has full SU(2) spin sym- 
metry but only a discrete charge symmetry. This dis- 
tinction leads, for instance, to the critical hybridization 
width Fci having a sublinear dependence on UcS [Eq. 
(j34p ] whereas its counterpart Fc2 is superlinear in Ucs 
[Eq. ((37)) ]. The equivalence of the critical points Cs 
and Cc under spin-charge interchange signals the emer- 
gence of a higher SU(2) isospin symmetry at both these 
renormalization-group fixed points. 

For r- > i, we find that (just as in the pseudogap 
Anderson modeP^) the symmetric strong-coupling fixed 
point of the pseudogap Anderson-Holstein model is un- 
stable with respect to any breaking of degeneracy be- 
tween the four impurity levels, i.e., to any UeS 7^ 0. As 
a result, the UeS-^ plane is divided into just two phases: 
local-moment for all t/off > and local-charge for all 
C/off < 0. 
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FIG. 10: Schematic renormalization-group flows on the Ueff- 
r plane for the particle-hole-symmetric pseudogap Anderson- 
Holstein model with a band exponent < r < i. Trajec- 
tories with arrows represent the flow of the couplings with 
decreasing temperature. Dashed lines connecting unstable 
fixed points (open circles) separate the basins of attraction of 
stable fixed points (filled circles). The right and left dashed 
lines represent the boundary values Fci and defined in 
Eqs. (|34p and (|37|l . respectively. The asymmetry of the fiows 
about the line UcS = stems from differing symmetries of 
the model in the spin and charge sectors. See the text for a 
discussion of each fixed point. 



IV. RESULTS: GENERAL MODEL WITH BAND 
EXPONENT < r < 1 

This section treats the pseudogap Anderson-Holstein 
model vifith a band exponent < r < 1 when either 
(the discrete) particle-hole symmetry is broken by a value 
6(1 = ed+\U 7^ or (the continuous) spin-rotation invari- 
ance is removed by a nonvanishing local magnetic field 
h [defined in Eq. It is found that increasing \5j\ or 
\h\ can drive the system from the local-moment or local- 
charge phase into one of several strong-coupling phases 
that are not present in the baseline case 5d = h = Q. 
The transitions between these phases take place at inter- 
acting quantum critical points in the same universality 
class as the asymmetric critical points of the pseudogap 
Anderson model. 

All numerical results presented in this section were ob- 
tained for an impurity with U = 0.5, for a bosonic energy 
Wo = 0.1, and for NRG discretization parameter A = 3. 

A. Phase boundaries 

1. Weak bosonic coupling 

Figure [11] plots phase boundaries of the pseudogap 
Anderson-Holstein model on the T-6d plane for zero mag- 
netic field, for band exponents r — 0.4 (left) and 0.6 
(right), and for three bosonic couplings A = 0, 0.1, and 
0.1414 that can all be associated [via Eq. ^] with ef- 
fective Coulomb interactions U > 0. A value 6d 
breaks particle-hole symmetry but leaves in place the 
SU(2) spin symmetry. The system remains in the local- 
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FIG. 11: (Color online) Phase boundaries of the zero-field 
pseudogap Anderson-Holstein model on the V-Sd plane for 
weak bosonic couplings A < Ao = 0.15812(1) and band ex- 
ponent r = 0.4 (left), r = 0.6 (right). Data calculated for 
U = 0.5, h — 0, ujo = 0.1, and the three values of A listed in 
the legend. 

moment phase for \6d\ < \tj and F < Fci(r, [/, 5^, A) = 
Fci(r, [/, —5d: A). Otherwise it lies in one of the strong- 
coupling phases described in Sec. Ill B II symmetric 
strong-coupling for 5d = 0, ASC_ for 5d > 0, or ASC+ 
for Sd < 0- Just as in the pseudogap Anderson modelfi^ 
the symmetric strong-coupling phase can be reached only 
for r < i; for r > ^, the symmetric strong-coupling 
fixed point is unstable and Fci(r, [/, (5^) diverges as Sd 
approaches zero. The contraction of the local-moment 
phase (i.e., the reduction of Fci) with increasing A and/or 
\Sd\ can be attributed to one or both of the Ud ^ 1 im- 
purity levels being drawn down in energy closer to the 
Ud = I ground states. This energy shift enhances in- 
direct spin-flip scattering between the ground states via 
the excited states and favors conduction-band quenching 
of the impurity degrees of freedom. 



2. Strong bosonic coupling 

Figure [12] plots phase boundaries of the particle-hole- 
symmetric pseudogap Anderson-Holstein model on the 
T-h plane for r — 0.4 (left) and 0.6 (right), and for 
three strong bosonic couplings in the range A > Aq. 
Here, Eq. ^ gives U < 0, and for local magnetic field 
h = the lowest-energy impurity states are spinless but 
have a charge Q = ±1. Application of a field h 
destroys SU(2) spin symmetry but preserves particle- 
hole symmetry. The field pulls one or other of the 
Q = excited states down in energy, thereby enhanc- 
ing virtual scattering between the Q = ±1 states. The 
model remains in the local-charge phase for \h\ < \U\ 
and F < Fc2(r, C/, /i. A) = Tc2ir,U,-h, \). Otherwise it 
lies in the symmetric strong-coupling phase (for r < ^ 
and h = 0) OT in one of two new asymmetric strong- 
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FIG. 12: (Color online) Phase boundaries of the particle-hole- 
symmetric pseudogap Anderson-Holstein model on the T-h 
plane for strong bosonic couplings A > Ao = 0.15812(1) and 
band exponent r — 0.4 (left), r = 0.6 (right). Data calculated 
for U = —2ed = 0.5, oJo = 0.1, and the three values of A listed 
in the legend. 



coupling phases: ASCj, with ground-state spin z compo- 
nent 5"^ = -i for /i > 0, or ASC^ with = ^ for h < 0. 
For 7" > i, Tc2{r, U, h, A) diverges as h approaches zero, 
a consequence of the instability of the symmetric strong- 
coupling fixed point in this range of band exponents. 

Although the two panels of Fig. [T^] look quite similar to 
their counterparts in Fig. 1111 there are small departures 
such as the absence from Fd versus 6d for r — 0.6 of 
points of inflection corresponding to those seen in F2c 
versus h. These differences presumably arise from the 
differing symmetries of the two cases. 



B. Impurity thermodynamic properties 



Weak bosonic coupling 



Figure [T^] plots the temperature dependence of 
TxsAmp, i^Xcimp, and Simp for r ^ 0.4, U = 0.5, 
Ed — —0.2 (or Sd — 0.05 > 0), A = 0.05, and seven values 
of F straddling Fd . The low-temperature limiting behav- 
iors in the local- moment phase (F < Fd) are identical to 
those found at particle-hole symmetry (see Fig. [5]). In 
the ASC_ phase (F > Fd), however, the T = proper- 
ties Txs,imp = liT'Xc.imp = -Simp = indicate complete 
quenching of the impurity degrees of freedom, in contrast 
to the partial quenching found in the symmetric strong- 
coupling phase. Exactly at F = Fd , the low-temperature 



properties Tx. 



s,imp 



0.1419, iTxc 



0.0147, and 



S'imp — 0.705 can be taken to characterize the critical 
point C_ separating the two stable phases. These prop- 
erties coincide with those found for (5d < at the critical 
point C-i- between the local-moment and ASC+ phases, 
and also with those of the asymmetric critical points of 
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FIG. 13: (Color online) Thermodynamic properties of 
the particle-hole-asymmetric pseudogap Anderson-Holstein 
model in zero magnetic field near the spin-sector critical point 
C-: Temperature dependence of the impurity contribution to 
(a) the static spin susceptibility Xs,imp multiplied by temper- 
ature, (b) the static charge susceptibility Xc.imp multiplied by 
temperature, and (c) the entropy S'lmp, for r = 0.4, Sd = —0.2, 
U = 0.5, ujo = 0.1, A = 0.05 < Ao ~ 0.158, and the seven val- 
ues of F — Fci labeled in the legend. Filled (open) symbols 
connected by guiding lines represent data in the local-moment 
(ASC_) phase, while thick lines without symbols show the 
critical properties at C_. Ns = 3 000 states were retained 
after each NRG iteration. 



the pseudogap Kondo or Anderson models;^* however, 
for r > r* ~ 3/8 these properties differ from from those 
of the corresponding symmetric critical point Cs. (Our 
observation that the asymmetric critical value of 
slig htly greater than In 2 is consistent with Refs. and 

EI) 

For F close to Fd , the crossover of the thermodynamic 
properties away from their critical values can be used to 
define a crossover scale that obeys Eq. ([55]) . Table HIl 
lists the correlation-length exponent i^i (r) obtained for 
five values of r. 



2. Strong bosonic coupling 



Figure [14] plots the temperature dependence of 
TXs.imp, I^Txcimp, and S'imp for r = 0.4, U = -2ed = 0.5 
(i.e., Sd = 0), h = 0.05, A = 0.2, and various F straddling 
the critical value Tc2- Just as was found at particle- hole 
symmetry, the thermodynamic properties for A > Aq are 
related to those for A < Aq by interchange of spin and 
charge degrees of freedom. For the cases shown in Fig. 
[TU it is the charge susceptibility that most clearly distin- 
guishes the critical point {jTxc.imp — 0.1419) from 
the local-charge and ASC^ phases {jTxcAmp = \ and 0, 
respectively). 
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FIG. 14; (Color online) Thermodynamic properties of the 
particle-hole-symmetric pseudogap Anderson-Holstein model 
in nonzero local magnetic field near the charge-sector critical 
point Cj,: Temperature dependence of the impurity contribu- 
tion to (a) the static spin susceptibility Xs.imp multiplied by 
temperature, (b) the static charge susceptibility Txc, imp mul- 
tiplied by temperature, and (c) the entropy Simp, for r — 0.4, 
U = -2ed = 0.5, h = 0.05, ujq = 0.1, A = 0.2 > Ao ~ 0.158, 
and the seven values of F — labeled in the legend. Filled 
(open) symbols connected by guiding lines represent data in 
the local-charge (ASC4.) phase, while thick lines without sym- 
bols show the critical properties at C^. A^s = 3 000 states were 
retained after each NRG iteration. 



C. Local response and universality class 

In the vicinity of the quantum critical points C± sep- 
arating the local-moment and ASC± phases, the local 
spin responses exhibit power-law behaviors described by 
Eqs. psp and (|42p . Table HIl lists critical exponents at 
C± obtained for five different values of r. To within their 
estimated accuracy, the exponents obey the hyperscaling 
relations Eq. (|43p , providing evidence for the interacting 
character of the critical points. 

Comparison between Tables U and |lT] shows that the 
symmetric critical point Cs and its asymmetric coun- 
terparts C± have the same low-temperature physics for 
r = 0.2 and 0.3, but not for r = 0.4. This pattern is con- 
sistent with the pseudogap Kondo model, where the Cs 
and C± critical points are identical for < ?- < r* ~ 3/8 
but distinct for r* < r < 1. In the latter range, the ex- 
ponents listed in Table HIl coincide to within small errors 
with those for the particle-hole-asymmetric pseudogap 
Kondo model given in Table II of Ref. [23l 

In the vicinity of the critical points C-|-4 marking the 
transitions from the local-charge phase to the ASC-f-.i 
phases, the local charge responses exhibit power-law be- 
haviors described by critical exponents equal (within 
small errors) to the local-spin exponents of the C± critical 
points of the pseudogap Kondo and pseudogap Anderson- 



TABLE 11: Exponents describing the local spin response at 
the critical points C± of the particle-hole-asymmetric pseu- 
dogap Anderson-Holstein model, evaluated for five band ex- 
ponents r. The critical exponents are defined in Eqs. psp 
and (|42p . A number in parentheses indicates the estimated 
random error in the last digit of each exponent. 



r 


1^1 


/3i 


l/Si 


Xi 


71 


0.2 


6.22(1) 


0.15(1) 


0.02630(2) 


0.9488(2) 


5.85(6) 


0.3 


5.14(1) 


0.34(1) 


0.07364(1) 


0.8629(3) 


4.41(3) 


0.4 


4.29(1) 


0.59(1) 


0.1569(1) 


0.7275(3) 


3.12(2) 


0.6 


1.78(1) 


0.188(1) 


0.1173(2) 


0.7896(4) 


1.41(1) 


0.8 


1.27(1) 


0.079(1) 


0.0644(5) 


0.879(1) 


1.10(1) 



Holstein models. We are led to conclude that these criti- 
cal points all lie in the same universality class. Given the 
different symmetries of the Anderson-Holstein model un- 
der spin and isospin rotation, the equivalence of the C± 
and C^,^ critical points under spin-charge interchange is 
a nontrivial finding, distinct from the equivalence of the 
Cs and Cc critical points at particle-hole symmetry. 



RESULTS: DOUBLE QUANTUM DOTS 
WITH [/a = 



In the pseudogap Kondo and Anderson models, the 
interacting quantum critical points found for band expo- 
nents < r < 1 are replaced for r > 1 by first-order 
QPTs that arise from renormalized level crossings be- 
tween spin-doublet and spin-singlet ground states of the 
impurity."2-5 Similar behavior is expected in the pseudo- 
gap Anderson-Holstein model. This section focuses on 
the particular case r = 2 that is of particular interest 
because it has a possible realization in double quantum 
dots. Below we present results not only for the impurity 
contributions to thermodynamic properties but also for 
the linear conductance of such a double-dot system in the 
vicinity of its spin- and charge-sector QPTs. 



A. Pseudogapped effective model for double 
quantum dots 

The motivation for focusing on the case r = 2 comes 
from theoretical studiea^Sr— of two lateral quantum dots 
coupled in parallel to left (L) and right (i?) leads, and 
gated in such a manner that the low-energy physics 
is dominated by just one single-particle state on each 
dot. It is assumed that one of the dots (dot 1) is small 
and hence strongly interacting, while the other (dot 2) 
is larger, has a negligible charging energy, and can be 
approximated as a noninteracting resonant level. This 
setup can be described by the two-impurity Anderson 
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Hamiltonian 

HdD = "J"'^ + Ulfll^fli^ + ££k4k<TC£ka 



(49) 



Here, dja annihilates an electron of spin z component a 
and energy Sj in the dot j {j = 1, 2), fij^ = dj^dja is the 
number operator for such electrons, and c^ko- annihilates 
an electron of spin z component a and energy Sik in 
lead £ (£ = L, R). For simplicity, the leads are assumed 
to have the same dispersion eik = £k corresponding to 
a "top-hat" density of states p{e) — pqQ{D — \e\) with 
Po = (^i*)^^, and to hybridize symmetrically with the 
dots so that VjL — VjR. Under these conditions, the 
dots couple only to the symmetric combination of lead 
electrons annihilated by Cko- — (cLku + CR]ia)/V2 with 
effective hybridization matrix elements Vj = a/2Vj£. 

A key feature of Eq. (|49p is the vanishing of the dot-2 
Coulomb interaction U2 associated with a Hamiltonian 
term C/2"-2t"-24.- This allows one to integrate out dot 2 to 
yield an effective Anderson model for a single impurity 
characterized by a level energy ei, an on-site interaction 
Ui , and a hybridization function^S 
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FIG. 15: (Color online) Phase diagrams of a f/2 = double- 
quantum-dot device with Ui = 0.5, Fi = 0.05, £2 = 0, F2 = 
0.02, and cjq = 0.1: (a) The A^-£i plane for h = Q. (b) The A^- 
h plane for £1 = —^Ui. Note the near-linearity of the phase 
boundaries in each panel when plotted against the square of 
the bosonic coupling. 



B. Phase boundaries 



ri(^) = Jir2 ^(D ~ kl), (50) 

[e ~ 62) +I2 

where Tj = irpoV^ for j — 1, 2. The presence of dot 2 in 
the original model manifests itself here as a Lorentzian 
hole in ri(e) of width F2 centered on £ = £2. For £2 = 
(a condition that might be achieved in practice by tun- 
ing a plunger gate voltage on dot 2), ri(£) cx £^ in the 
vicinity of the Fermi energy, providing a realization of 
the r — 2 pseudogap Anderson modeli^ 

In the remainder of this section, we consider the 
double-dot device introduced in Ref. 20l, augmented by a 
Holstein coupling between dot 1 and local bosons. Such 
a system, modeled by a Hamiltonian Hdd + wgo^a -I- 
A(ni — l)(a -I- a^), can be mapped (following Ref. [20I ) 
onto the effective single-impurity model 

O" k,(T 

+ ^ (4ck, + H.c.) + A(fii - 1) (a + at) (51) 

k,cr 

with the hybridization function ri(£) = 
ttN^'^ Ek ViH>^ - Ek) as defined in Eq. §0^. 

All numerical results presented in the remainder of this 
section were obtained using the effective one-impurity 
pseudogap Anderson-Holstein model [Eq. (|5ip ] for a 
strongly interacting dot 1 having Ui — 0.5 and for a 
bosonic frequency ujq =0.1. The NRG calculations were 
performed for a discretization parameter A = 2.5. 



1. Weak bosonic coupling 

Figure [TST a) shows the phase diagram of the U2 — 
double-quantum-dot device, as mapped to the pseudo- 
gap Anderson-Holstein model, on the A^-£i plane in the 
absence of any magnetic field. For A < Aq — 0.1582848, 
decreasing the dot-1 energy starting from a large pos- 
itive value drives the system from the ASC_ phase to 
the local-moment phase (LM) at £1 = ef^ < 0, and 
then from the local-moment phase to the ASC+ phase at 
£1 = = —Ui — £]^^. For A > Ao, the system instead 
lies in the ASC_ phase for all £1 > —^Ui, in the ASC+ 
phase for all £1 < —^Ui, and in the local-charge phase 
(LC) only along the line £1 = —Ui/2 of strict particle- 
hole symmetry. 

One of the most notable features of Fig. [TSl a) is the 
linear dependence of ef^ on A^, which implies a linear 
dependence on the polaron energy £p = X^/loq. Since 
£^j, = —Ui — ef_^, it suffices to focus on the phase bound- 
ary between the ASC_ and LM phases. In the atomic 
limit F = 0, one expects this boundary to be defined by 
the degeneracy of the n^j = and n^i = 1 impurity levels, 
i.e., the point where the renormalized dot-1 level energy 
[cf. Eq. ([TT]) ] satisfies £1 = £1 -I- £p = 0, a condition that 
implies e^^ — —Sp. 

The location of the phase boundary for F > can 
be estimated using the poor-man's scaling equations dis- 
cussed in Sec. Ill Cl Equation (pij) implies that for r > 1, 
the effective value of the dimensionless scattering width 
r/D = (D/ DY~^{T / D) decreases monotonically under 
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reduction of the half-bandwidth from D to D. If F < D, 
this decrease rules out the possibility of entry into the 
mixed-valence regime under the criteria laid out at the 
end of Sec. Ill CI Moreover, the decrease of F is so rapid 
that any entry to the local-moment regime and subse- 
quent mapping to the pseudogap Kondo problem will 
yield a sub-critical exchange coupling, placing the system 
in the local-moment phase. Under these circumstances, 
the boundary between the local-moment phase and the 
asymmetric strong-coupling phase ASC_ [see Sec. IIIB II 
and in particular Fig. [IJb)] is effectively determined by 
the condition ed{D — 0) — for a level crossing between 
the renormalized energies of the empty and singly occu- 
pied impurity configurations. With the approximation 
U = oo, and using the expansion 

S(a, x) ~ 1 — + O(x^) for x > 0, (52) 

a + 1 

Eq. (flT)) can be integrated to yield 

r 

4,c--^-"iep' (53) 

with 

/ ujo . D + ujo\ F 

The predicted value ai = 0.986 is in good agreement 
with the one ai — 0.988 that describes NRG results for 
the Anderson-Holstein model with U = 0.5, F — 0.05, 
and a pure r — 2 power-law hybridization function (data 
not shown) . The phase boundary for the mapped double- 
quantum-dot system plotted in Fig. [TSla) can be fitted 
with a reduced value ai = 0.80 that can be attributed 
to the fact that the hybridization function ri(e) in Eq. 
([SO]) assumes a power-law form only for |e| <C A2 ^ -D. 

2. Strong bosonic coupling 

Figure [rsl b) shows the phase diagram of the U2 — 
double-quantum-dot system on the X^-h plane at 
particle- hole symmetry. For A > Aq, decreasing the local 
magnetic field from a large positive value takes the sys- 
tem from the ASCj. phase to the local-charge phase (LC) 
a.t h — he > 0, and then from the local-charge phase to 
the ASC-|- phase at h — —he- For A < Aq, by contrast, 
the system is in the ASC4. or ASC-|- phase for /i > or 
h < 0, respectively, and in the local-charge phase only 
along the line h — 0. 

The phase boundaries shown in Fig. llSf b) are nearly 
linear in A^ — Aq or, equivalently, linear in 2ep — U. In 
the atomic limit F = 0, one expects the boundary be- 
tween the ASC4. and LC phases to be located at the 
point where the singly occupied Sz = —■^ impurity state 
crosses energies with the degenerate pair of impurity 
states having n^; = 0, 2, i.e., to satisfy — = ■^U or 
he — 2£p — U. For F > 0, a generalization of the poor- 
man's scaling analysis of Sec. Ill Ci to incorporate the local 
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FIG. 16: (Color online) Crossover temperature scales near 
QPTs in a (72 = double-quantum-dot device with Ui =0.5, 
Fi = 0.05, £2 = 0, r2 = 0.02, and cjo = 0.1: (a) vs 
Aei — £1 — ef^ in the ASC- phase for h = and A = 0.1. 
(b) T2 vs Ah — h~ he in the ASCj. phase for ei = —0.25 and 
A = 0.2. Straight lines represent power-law fits, which yield 
the correlation- length exponents = 1/2 = Ns = 1000 
states were retained after each NRG iteration. 



field would be expected to yield corrections to this result 
along lines similar to the corrections found in the regime 
A < Aq. Empirically, we find that the data in Fig. ITSTb) 
for < A — Ao ^ Ao can be fitted to the form 

he a2i2ep ~ U) (55) 

with a2 = 0.69. However, a larger value Q!2 — 1 is re- 
quired to describe the data points for A ~ 2Ao, indicating 
that the critical field is not strictly linear in 2£p — U. 



C. Crossover scales 

As discussed in Sec. IIII C| we can use the NRG 
spectrum to identify temperature scales characterizing 
crossovers between different renormalization-group fixed 
points. Figure [UK a) plots the crossover scale T* versus 
Aei = El — on approach to the local-moment phase 
boundary from the ASC_ phase for a weak bosonic cou- 
pling A = 0.1, while Fig. fTHT bi shows the scale T2 versus 
Ah = h — he in the ASC^ phase near the local-charge 
phase boundary for A = 0.2. Tj* and Tj* vanish at the 
phase boundaries in the manner 

T* \e - e+^l"' as e ^ 4„ (56) 

and 

T; (x\h- hel"^ as h^ he, (57) 

with correlation- length exponents vi —1^2 — 1- This 
linear vanishing of crossover scales is consistent with 
the level-crossing nature of the QPTs of the pseudogap 
Anderson-Holstein model for r = 2. 
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D. Impurity thermodynamic properties 

As in the cases r < 1 considered above, the temper- 
ature variation of the impurity contributions^ to the 
static spin and charge susceptibihties and to the entropy 
can be used to distinguish the strong-couphng phases of 
the U2 — double-quantum-dot system from the phases 
with residual local spin or charge degrees of freedom. 



1. Weak bosonic coupling 

Figure [17] plots the temperature dependence of the im- 
purity thermodynamic properties Txs.imp, jl^Xc.imp, and 
•S'imp for h — 0, a weak bosonic coupling A = 0.1, and 
various values of ei straddling the upper transition. For 
El = ~ —0.124985 (lines without symbols in Fig. 
flT)) . the low-temperature limiting values Txs, imp = 1/6, 
jTxc,imp = 1/18, and S'imp = ln3 are those character- 
istic of the valence- fluctuation fixed point: the point of 
degeneracy between impurity occupancies = and 
= I corresponding to Eq. (jSip with effective couplings 
El = Fi = A = and Ui = 00. If £1 deviates slightly from 
its critical value, the properties trace their critical behav- 
iors at high temperatures, but cross over below a scale T* 
to those either of the local-moment phase, where there 
is a residual spin-i degree of freedom (Txs,imp = 1/4, 
Txc,imp — 0, and S'imp = In 2) or of the ASC_ phase 

(-^Xs,imp — ^Xcjimp — S'imp — 0)- 



2. Strong bosonic coupling 

Figure [m shows Txs.imp, I^Xcimp, and Simp vs T at 
particle-hole symmetry (ei — —^Ui) for a strong bosonic 
coupling A = 0.2 and various local magnetic fields h 
straddling the critical value — 0.284959. Here, in con- 
trast to Fig. [T71 Txs,imp falls to zero in both the local- 
charge phase and the ASC^ phase, signaling the suppres- 
sion of spin fluctuations at the impurity site. However, 
the flows of Txcjimp with decreasing temperature clearly 
reveal the existence of a QPT separating the ASC^ and 
local-charge phases. Exactly at the critical value h — he 
(lines without symbols in Fig. [TS|) . jTxc.imp is pinned 
at low temperatures at the value 1/6 expected at the 
point of degeneracy between the empty, spin-down, and 
doubly occupied impurity configurations. For h deviat- 
ing slightly from he, Txc,imp traces the critical behavior 
at high temperatures but eventually crosses below a scale 
T2 to a limiting value of cither 1 in the local-charge phase 
or in the ASCj. phase. 

The temperature dependence of the spin (charge) sus- 
ceptibility in Fig. [TH] mirrors that of the charge (spin) 
susceptibility in Fig. [iTl By contrast, the behavior of 
Simp (7") is equivalent in the two cases. These properties 
suggest that, as found for the interacting quantum crit- 
ical points for band exponents < r < 1 (Sees. IIIII and 




FIG. 17: (Color online) Thermodynamic properties of a 
U2 ~ double-quantum-dot device near a spin-sector QPT: 
Temperature dependence of the impurity contribution to (a) 
the static spin susceptibility Xs.imp multiplied by temperature, 
(b) the static charge susceptibility Xc.imp multiplied by tem- 
perature, and (c) the entropy Simp, for Ui = 0.5, Fi = 0.05, 
£2 = 0, r2 = 0.02, ft = 0, wo = 0.1, A = 0.1, and nine val- 
ues of |ei — £^^| where ef^ ~ —0.124985. Properties at the 
transition (thick lines without symbols) are those expected at 
a level crossing between the local-moment phase (filled sym- 
bols) and the ASC_ phase. Ns — 3 000 states were retained 
after each NRG iteration. 



lIVp . the quantum phase transitions into/out of the LC 
phase at 6d = take place at points of enhanced symme- 
try where the system acquires an SU(2) isospin invari- 
ance to match the global SU(2) spin invariance of the 
Anderson-Holstein Hamiltonian in zero magnetic field. 

Both at weak and strong bosonic couplings, the fact 
that the impurity properties in the quantum-critical 
regime are those of the valence-fiuctuation fixed point (or 
its analog under interchange of spin and isospin) is en- 
tirely consistent with the picture of each QPT as arising 
from a renormalized level crossing. Moreover, crossover 
scales Tj* and Tj* extracted from the thermodynamic 
properties are identical up to a constant multiplicative 
factor to those identified from the NRG spectra (see Sec. 

EC]). 



E. Linear conductance 

It is generally impractical to measure the impurity 
thermodynamic properties of a quantum-dot device. 
Rather, the primary experimental probe of lateral quan- 
tum dots is electrical transport. The linear conductance 
of the boson-coupled double-quantum-dot system mod- 
eled by Eq. can be calculated from the Landauer 
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FIG. 18: (Color online) Thermodynamic properties of a 
U2 = double-quantum-dot device near a charge-sector QPT: 
Temperature dependence of the impurity contribution to (a) 
the static spin susceptibility Xs,imp multiplied by tempera- 
ture, (b) the static charge susceptibility Xc.imp multiplied by 
temperature, and (c) the entropy Simp, for Ui = —2ei = 0.5, 
Ti = 0.05, £2 = 0, Ta = 0.02, cjq = 0.1, A = 0.2, and nine 
values of \h — hc\ where he — 0.284959. Properties at the tran- 
sition (thick lines without symbols) are those expected at a 
level crossing between the local-charge phase (filled symbols) 
and the ASC4. phase. Ns = 3 000 states were retained after 
each NRG iteration. 
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FIG. 19: (Color online) Linear conductance of a. U2 = 
double-quantum-dot device near a spin-sector QPT: (a) Lin- 
ear conductance g vs Aei — £1 — ef^ for Ui — 0.5, Fi — 0.05, 
£2 = 0, F2 = 0.02, /i = 0, uo = 0.1, A = 0.1, and dif- 
ferent temperatures T specified in the legend as multiples of 
Tko = 7 X 10"". The retention of iV^ = 1 000 states after each 
NRG iteration accounts for the small discrepancy between 
e;^'^ ~ -0.1249871 and its value in the case Ns = 3 000 shown 
in Fig. 1171 (b) The same data scaled as {g ~ gc)TKo/T vs 
Aei/T, where gc{T) is the conductance at ei = e^^. The solid 
line was obtained from Eq. (|62|) by approximating Aq{x) = aq 
and using values of a, and ijjq/{£i ~ e^) fitted from Aii{lu). 



formula 

2 foo 

9iT) = T E / '^^ {-df/duj)[-lm%{^)], (58) 

where f{uj,T) ~ [exp(a;/T) + 1]^^ is the Fermi-Dirac 
distribution function and Xt(w) — npo Vi Gfj{uj) Vj 

with Gf^iio) = -*/o°°dte-*({d,,(t),4(0)}). For U2 = 
0, one can re-express^i 

-ImTaiuj) = [1 - 27rF2/92(w)]7rr(cj)A!^i(w) +7rF2P2(w) 
+ 2n{uj-e2)r{uj)p2{uj)ReG^,{cj), (59) 

where T{uj) is as defined in Eq. ([SHI) . P2(e) — (r2/7r)[(£ — 
£2)^ + r^]"^ is a Lorentzian of width F2 centered on en- 
ergy £2, and Aii{uj) = —Tr~^lmGii{u}). All quantities 
entering Eq. ((58|) are known exactly with the sole excep- 
tion of Gii{u!), the full dot-1 spin-cr local Green's func- 
tion Gii{uj) taking into account both electron-electron 
{Ui) and electron-boson (A) coupling. 

We have used standard NRG methods^^ to obtain the 
dot-1 spectral function AJ]^(w) from the effective one- 
impurity pscudogap Anderson-Holstein model described 
by Eq. ([^Tjl . After KeGii{u!) has been obtained via 
the Kramers-Kronig relations, Eqs. ((58)) and (|59)) yield 
the linear conductance. We present results only for zero 
magnetic field [where A|]^(lli) = Af]^(a;)] and/or for strict 



particle-hole symmetry [where ^|;^(a;) = ^{;^(— w)], spe- 
cial cases in which the up and down spin channels con- 
tribute equally to the conductance. Temperatures are 
expressed as multiples of Tko = 7 x 10"'', the Kondo 
temperature for the conventional (i.e., metallic or r — 0) 
one- impurity Anderson model with U = —2ed = 0.5 and 
F — 0.05, which serves as a characteristic scale for the 
many-body physics of the problem. 

1. Weak bosonic coupling 

Figure fTOT a) plots g versus A£i — ei ~ e^^^ for a weak 
bosonic coupUng A = 0.1 and five temperatures T Hsted 
in the legend. At T = 0, the linear conductance g is struc- 
tureless and takes its maximum possible value 2e^/h, sig- 
naUng perfect electron transaction through the system. 
However, at T > 0, g versus A£i develops clear min- 
ima on either side of a maximum located precisely on 
the boundary A£i = between the local-moment and 
ASG_ phases. The peak-and-valley structure becomes 
more prominent upon increasing temperature up to sev- 
eral times Tko, making it amenable to experimental ob- 
servation. Similar features have been reported [see Fig. 
2(b) of Ref . [2l| for a double-quantum-dot system without 
bosonic coupling. 

The essential features of the results shown in Fig.fTWa) 
can be understood from the fact that — ^just as in the 
case of zero bosonic coupling^! — near the QPT, the low- 
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energy part of the dot-1 spectral function is dominated 
at low temperatures by a quasiparticle peak at frequency 
iijq (X El — El . Upon raising the temperature, this peak 



rapidly disappears once T > uiq 
behavior by 



We approximate this 



^?i(cj) ~ v4,K/r) 6{lj - uj,) for « T^, (60) 

where Aq (x) is an unknown scaling function that satisfies 
Ag{x) for x < 1 and Aq{x) a, > for x > 1. 
Hilbert transformation of Eq. (|60)) leads to 

Re Gl, iu,T)c^Ro{ujq,T) + Riiujq,T)^ + M^IIl , 

J- 2 LUq — OJ 

(61) 

where Rq and i?i are determined by the form of A^i (w) at 
I'^l ^ ^2, and may vary with Si and hence uig. Inserting 
Eqs. dini) and dM]) into Eqs. ^ and §^ yields, for 
£2 = andT<r2, 



2e 



9 = 



h {^^^[^2 

Aq{Uq/T)^-^ 



—— I dx 



a;2 - {ujq/TY {e^ + 1)2 



O ^ 



(62) 

4, 



The first line in Eq. ((62|) . which describes resonant tun- 
neling through dot 2, dominates the conductance both 
for T <^ ujq and for T 3> Wg. However, for T ~ Wq, the 
conductance is dominated by the first term in the square 
brackets, which arises from the A^i term in Eq. (|59| . 
To good approximation, the conductance near the QPT 
(where ujq = 0) can be expressed as 



g{ei,T)=gie+T)+Tg, 



£l 



T 



(63) 



with (7i(0) — 0. Figure [TWb) shows that this form is 
obeyed well by the NRG results. The precise scaling 
function gi cannot be determined without knowledge of 
Aq{x)^ but the zeroth-order approximation Aq(x) — aq 
produces a reasonably good description of the numeri- 
cal data. This scaling collapse of the finite-temperature 
conductance feature provides a clear signature of the un- 
derlying T — Q phase transition that may be sought in 
experiments. 



2. Strong bosonic coupling 

As one would expect given the equivalence under spin- 
charge interchange of the thermodynamic properties at 
the LM-ASC± and LC-ASC-f-^ phase boundaries, the 
variation of the conductance with h around the criti- 
cal field he for A > Aq is very similar to the varia- 
tion of g with £i near e^^,. The system exhibits perfect 



electron transmission {g = 2e'^/h) at T = and with 
increasing temperature develops an increasingly promi- 
nent peak-and-valley signature of the QPT. This signa- 
ture can be understood as arising from the existence of 
quasiparticle peaks A|j^(w) ~ Aq{ujq/T) 6{u! — ojq) and 
yl][^-^(a;) ~ Aq{u!q/T) 6{u!+ujq) at a frequency cuq cx h — hc- 
Analysis similar to that applied in the case of weak 
bosonic couplings leads to the prediction 



g{h,T)^g{h,,T) + T~g^ 



T 



(64) 



a scaling that is indeed displayed by the numerical data. 
(We do not show these data explicitly due to their simi- 
larity with Fig.IHl) 



VI. SUMMARY 

We have conducted a study of the pseudogap 
Anderson- Holstein model describing a magnetic impurity 
level that hybridizes with a pseudogapped fermionic host 
with a density of states vanishing as |e|'" at the Fermi en- 
ergy (e = 0), and that is also coupled, via its charge, to 
a local-boson mode. The reduction of the density of low- 
energy band excitations leads to quantum phase transi- 
tions (QPTs) that can be classified into different types 
depending on the strength of the impurity-boson cou- 
pling and on the presence or absence of particle- hole and 
time-reversal symmetry. The main results are as follows: 

(1) Under conditions of strict particle- hole and time- 
reversal symmetry, the pseudogap Anderson-Holstein 
model with exponent < r < ^ features two types of 
continuous QPT. For a weak (strong) impurity-boson 
coupling that results in a positive (negative) effective 
Coulomb interaction between electrons in the impurity 
level, increasing the impurity-band hybridization from 
zero drives the system through a continuous QPT be- 
tween a local-moment (local-charge) phase, in which a 
two-fold degree of freedom survives to T = 0, and a 
symmetric strong-coupling phase in which the impurity 
degree of freedom is quenched by the conduction band. 
Critical exponents characterizing the response to a local 
symmetry-breaking field suggest that these QPTs belong 
to the same universality class as the QPT of the particle- 
hole-symmetric pseudogap Anderson model. 

(2) For r > i, the symmetric strong-coupling fixed 
point is unstable (just as in the pseudogap Anderson 
model without bosons) and for weak (strong) impurity- 
boson couplings, a system exhibiting particle-hole and 
time-reversal symmetry always lies in the local-moment 
(local-charge) phase. 

(3) For weak impurity-boson couplings and away from 
particle-hole symmetry, the symmetric strong-coupling 
phase is replaced by two asymmetric strong-coupling 
phases, one corresponding to an empty impurity level 
and the other to double occupation of the impurity 
site. These phases are separated from the local-moment 
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phase by QPTs in the same universality class as those 
of the particle-hole-asymmetric pseudogap Anderson 
model. These QPTs are continuous and interacting for 
< r < 1, but first order for r > 1. 

(4) For strong impurity-boson couplings and in the 
presence of a magnetic field, the local-charge phase is 
separated by QPTs (again in the asymmetric pseudogap- 
Anderson universality class, and continuous for r < I 
but first-order for r > 1) from two asymmetric strong- 
coupling phases corresponding to single occupation of the 
impurity level with either a spin-up or a spin-down elec- 
tron. 

(5) For r = 2, the pseudogap Anderson-Holstein model 
provides a description of two quantum dots connected 
in parallel to current leads, where one dot is tuned to 
lie in a Coulomb blockade valley and is coupled via its 
charge to a local-boson mode, while the other dot is tuned 
to be effectively noninteracting and in resonance with 
the leads. The setup exhibits voltage- or magnetic-field- 
tuned QPTs of the level-crossing type. These QPTs pro- 
duce peak-and-valley features in the linear conductance 
that become more prominent upon increase of the tem- 
perature. Moreover, in the vicinity of the transitions, 
the conductance data collapse to a single function of the 
ratio of a symmetry breaking field to the absolute tem- 
perature. 
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Appendix A: Derivation of Poor-Man's Scaling 
Equations 

In this appendix, we outline the derivation of the poor- 
man's scaling equations p^ - (fT8l) discussed in Sec. Ill CI 
For this purpose, it proves convenient to work with the 
Anderson-Holstein Hamiltonian in the form 



where the hybridization functions 



~^ -^imp-band + -^imp-boson , 

where i?band, -ffboson, and -ffimp-boson are as defined in 
Eqs. (pb| . (Pc)) . and respectively, but Hiaip in Eq. 
Zal) is rewritten in more conventional fashion as 



(A2) 



and -ffimp-band in Eq. (|2dp is generalized to 



TTl 

imp-band 



k,(T 



r.(£) 



—y 



Vr^l^Sis - e^) = Tr\s/D\- eiD ~ \e 



(A4) 

for T = 0, 2 have the same power- law dependence as r(e) 
defined in Eq. Q. At the bare Hamiltonian level, one 
expects the hybridization matrix element Vb,k between 
the empty and singly occupied impurity configurations to 
be identical to that V2,k between the singly occupied and 
doubly occupied impurity configurations. However, this 
degeneracy can be broken under the scaling procedure. 

A canonical transformation H' ^ H' = e^H'e"^ with 
S as defined in Eq. (O yields 



H' = Hi 



imp 



^^band + -ffboson + ^imp-band ' (^5) 

where H[^p contains shifted parameters U [Eq. ^] and 
Ed [Eq. dll])], and 



k.cr 



V2Mnd^-a]dlck^ + U.c.}, (A6) 



H' 

imp-band 



with B as defined in Eq. ([iT 

We analyze the problem using a basis of many-body 
states composed as direct products of (i) fermionic states 
formed by the action of creation and annihilation opera- 
tors on \FS), the half- filled Fermi sea having Nk elec- 
trons of energy £k < 0, and (ii) occupation number 
eigenstates \n) of the transformed boson mode defined 
in Eq. ([14]). Since real occupation of states |n) with 
71 > is negligible in the anti-adibatic regime, we fo- 
cus on the states |0,0) = \FS) ^ |0), |cr,0) = 4|0,0), 
and 1 2,0) = a dl\ — a,0). Neglecting for the moment the 
effect of the hybridization [^^i'mp.band Eq. (|A6p ]. the 
energies of these states are denoted Eq, Ei — Eq + Ed, 
and E2=Ei+ed + U = 2Ei- Eo + U. 

We now consider the effect of an infinitesimal reduc- 
tion in the half-bandwidth from D to D ^ D + dD, where 
dD < 0. The goal is to write a new Hamiltonian H' sim- 
ilar in form to H' but retaining only conduction-band 
degrees of freedom having energies |£k| < D and having 
parameters id, U, and r,- adjusted to account perturba- 
tively for the band-edge states that have been eliminated. 

Let be the set of wave vectors k describing particle- 
like states having energies D < < D , and be the 
set of wave vectors describing hole-like state with energies 
—D < £k < —D. Tunneling of an electron from a K~ 
state into the empty impurity level, accompanied by the 
creation of n = 0, 1, . . . local bosons, transforms the state 
|0,0) to 



|0,0) = |0,0) 



H.C.}, 



(A3) 



OO ^ 



kGK-,o- 

(A/o^o)" 



n=0 



y/n\ jekj +Ei- Eo + nujQ 



c,,,\a,n)+0{V') 
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with energy 
Eo = Eq — 



N, 



oo ^ 

E- 

Z — / r) 



(A/u;o)2" 



n=0 



n! |ek| + El - Eq + noja 



+ 0{V'^). (A8) 



Here, 0(F™) schematically represents all processes in- 
volving at least m factors V^n.ki • • • ^^r„,k„- The deriva- 
tion of Eqs. (|A7I) and (lASp makes use of 



(n|e±"('''-'')|0) (n|e±"(^'-'')|0) 
(n|e^"" e 



= P-"V2 rr7lp±"5t^=Fa6|Q^) 



V2 



(±a)". 



(A9) 



Since N^^ Ekeif± K% - 7r-ir^(±L') S{ewTD), one can 
re-express the perturbed energy 



Eo-Eo- \dD\ 



2Toi-D) 



0{V% 



(AlO) 



where £{E) is the energy function defined in Eq. (fT9|) . 

Similarly, tunneling of an electron from the doubly oc- 
cupied impurity level into a K'^ state transforms |2, 0) 
to 



|2,0) = |2,0) 



E 



2,k 



oo ^ 



keA'+,(T 

(A/o.o)" 



Vn! Ek - - £d + nujo 



(All) 
0(^2) 



with energy 



Eo — Eo 



oo -j^ 

E II 



n=0 



\2 / 2 

-A /lu„ 

keK+,0- 
n! e\^ — U — ed + nuiQ 



Eo 



\dD\ 



2r2p) 



^£{D~U~ed) 



0(^3 



(A12) 



Finally, tunneling of an electron into the singly occupied 
impurity from a K~ state or from the singly occupied 
level into a state transforms |cr, 0) to 



k,0) = |a,0)- 



E 



2,k 



F — ■ 



L keAT- 

(A/c^o)" 



y/n\ kkl + i?2 - + ?^wo 



Ck, 



|2,n) 



(A/o.o) 



0,n) 



with energy 



-El — El 



E 

keA- 

(AM)'" 



2,k| 



n! |e 

' E 

k6A'+ 

: £;i - \dD\ 
+ 0{V^) 



v.\ + E2 — El 
00 _^ 

n=0 



(A/^o)'" 



So 



r2(-7?) 



ro(^) 



(A14) 



The 0(T/^) terms in each of the above states |(/), 0) in- 
clude terms to enforce normalization, i.e., (0, 0|(/>, 0) = 
(0,010,0) = 1. 

The modified energies can be used to define effective 
Hamiltonian parameters id = Ei — Eq and U ^ E2 + 
Eq — 2Ei. At the same time, for each k in the retained 
portion of the band (i.e., satisfying |ek| < D), Vo,k must 
be replaced by 



Vqu = 



/iVfc (0,0|Btcka-ffV,0) forek>0 
-VN^ (MIB^^TJ'IOTO) for ek < 0, 



(A15) 



and V2M must be replaced by 

-rrViVfc (<^\B^Ck,^^H'M for £k > 
•yiVfe (270|b4 _^7?'M) for ew < 0. 



l/2,k- 



It is straightforward to show that 

VrM = K,k + o(y3). 



(A16) 



(A17) 



We shall not attempt to evaluate the leading corrections, 
which turn out to be negligible in pseudogap (r > 0) 
cases. 

The infinitesimal band-edge reduction described in the 
previous paragraphs can be carried out repeatedly to re- 
duce the half-bandwidth by a finite amount from D to 
D < D. Equations (|A10|) and (|A14[) indicate that during 
this process, the impurity level energy evolves according 
to the scaling equation 



did 

db 



2Tq 



£{D-id) £{D + id) £{D + U+id) 

(A18) 

where Tt^± is the value of the rescaled hybridization func- 
tion at the reduced band edges e = ±D. Taking into 
account Eq. (jA12p as well, one sees that the on-site re- 
pulsion follows 



dU 
dD 



- 0,+ 



r, 



0: + 



2,+ 



£{D+id) £{D~id 



£{D + U+id) 



£{D~U^id) 
0(y3). (A19) 
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The band-edge hybridization functions rescale both 
due to the replacement of D by l) in Eq. ^ and due to 
the perturbative corrections to VV,k in Eq. ()A17| , leading 
to the scaling equation 

^£^=r^ + 0{V*). (A20) 
dD D 

The bare hybridization functions specified in Eq. (jA4|) 



are such that Tt,±{D) = T. For r > 0, moreover, Eq. 
(|A20I) shows that the band-edge hybridization functions 
are irrelevant (in the RG sense), and so too must be any 
differences among the renormalized values of the four hy- 
bridization widths. It is therefore an excellent approxi- 
mation to set Tr,± = F from the outset, leading to the 
simplified scaling equations given in Eqs. (I16l) - (|18p . 
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